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The  objectives  of  this  study  were  several  fold:  1)  To  examine  space  and  time 
dependent  transport  In  gallium  arsenide,  aluminum  gallium  arsenide,  and  Indium 
gallium  arsenide;  2)  To  begin  a  reformulation  of  quantum  transport  for  mixed 
states.  The  technique  used  for  studying  transient  transport  was  through 
solutions  to  the  moments  of  the  Boltzmann  transport  equation.  The  quantum 
transport  was  primarily  formulatory. 

In  examining  the  binary  and  ternary  transport  application  of  Vegards  rule  was 
applied  to  determining  the  effective  masses  and  a  variety  of  other  band  struc¬ 
ture  parameters.  The  parameters  for  Indium  arsenide  and  aluminum  arsenide  were 
obtained  from  previous  calculations  of  Littlejohn,  et.al.,  and  modified  for 
implementation  using  the  moments  of  the  Boltzmann  transport  equation.  The 
initial  calculations  for  steady  state  were  performed  using  parabolic  bands.  It 
became  Immediately  apparent  that  the  parabolic  band  structure  would  be  Inade¬ 
quate  for  examining  transport.  Thus  during  the  Initial  phases  of  the  study, 
the  moments  of  the  Boltzmann  transport  equation  were  modified  to  Include  the 
presence  of  generic  energy  versus  momentum  dispersion  relationships.  All 
calculations,  however,  were  performed  for  the  situation  In  which  e  +  > 

-h^k^/Zm*,  where  Cg  is  the  separation  between  the  conduction  band  mlnlmla 
and  the  top  of  the  valence  band.  Nonparabollclty  was  used  for  the  T,  L,  and  X 
portions  of  the  conduction  band  for  all  semiconductors  of  Interest. 

In  all  calculations  a  displaced  Maxwellian  distribution  function  was  used. 

For  the  case  of  parabolic  bands,  the  derived  moments  of  the  Boltzmann  transport 
equation  are  direct  and  have  been  published  In  variety  of  places  Including 
those  by  the  Principal  Investigator.  A  new  feature  Introduced  under  the  ARO 
study  was  the  Implementation  of  the  displaced  Maxwellian  for  nonparabollc 
bands.  In  this  case,  it  was  necessary  to  expand  the  formal  displaced 

Maxwellian,  f(r,it,t)  In  powers  of  k.  The  expansion  as  performed  In  this  study 
consisted  of  four  terms: 


f(r,ic,t)  =  fo  +  fi  +  f2  +  ^3 


(1) 


higher  order  terms  were  Ignored.  Most,  If  not  all,  previous  discussions  of 
nonparabollc  transport  have  Included  only  the  fg  and  fj  terms.  The  departure 
In  this  study  was  the  Inclusion  of  the  f2  and  f3  terms.  The  first  question 
investigated  was:  How  Important  is  the  Inclusion  of  the  f2  and  f3  terras? 
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To  determine  Che  Importance  of  the  £2  £3  terms,  Che  moments  of  Che 

Boltzmann  transport  equation,  which  had  already  been  obtained  exactly  for 
parabolic  bands  was  rederlved  for  parabolic  bands  using  the  four  term  expan¬ 
sion.  It  was  determined  Chat  the  form  of  the  moment  equations,  as  obtained 
using  the  complete  parabolic  displaced  Maxwellian,  and  Che  four  term  expansion, 
were  Identical  only  when  all  four  expansion  terms  were  Included.  This  result 
Indicated  Chat  all  previous  nonparabollc  studies  using  the  moment  equations  and 
only  the  fg  and  fj  terms  were  significantly  Incomplete.  The  summary  In  the 
next  few  paragraphs  discusses  the  significance  of  these  latter  contributions, 
as  well  as  the  general  nonparabollc  contributions  which  have  not  been  discussed 
In  the  past. 

To  begin.  It  Is  worthwhile  to  restate  the  obvious,  particularly,  when  nonpara¬ 
bollc  transport  Is  considered:  The  transient  dynamics  of  carriers  both  under 
uniform  and  nonuniform  field  conditions.  Involves  momentum  transport  rather 
than  velocity  transport.  The  moments  of  the  Boltzmann  transport  equation  de¬ 
scribe  the  time  dependence  of  carrier  momentum,  and  all  overshoot  phenomena 
which  Is  generally  ascribed  Co  velocity.  Is  rather  momentum  cmremhoot.  Under 
uniform  field  conditions,  the  time  rate  of  change  of  carriers  due  to  Inter¬ 
valley  transfer.  Is  affected  only  through  nonparabollc  contributions  to  the 
scattering  Intervals.  Under  nonuniform  field  conditions  the  divergence  term, 
dlv(nhk/m),  which  appears  In  the  continuity  equation.  Is  unchanged  In  form  for 
both  parabolic  and  nonparabollc  bands.  The  significance  physics,  however.  Is 
altered:  For  parabolic  bands  -ftk/m  Is  the  mean  carrier  velocity;  for  nonpara¬ 
bollc  bands  It  Is  not. 

With  regard  to  momentum  balance  for  uniform  fields,  the  form  of  Che  equation  Is 
the  same  for  both  parabolic  and  nonparabollc  bands.  The  only  alteration  Is  In 
the  nonparabollc  contribution  to  the  momentum  scattering  Integrals.  However, 
for  nonuniform  fields  there  Is  a  significant  alteration  In  the  pressure  Censor. 
The  alteration  In  the  pressure  tensor  for  nonparabollc  bands  arises  from  Che  f2 
contribution  In  the  expansion  of  the  distribution  function.  Complete  cancel¬ 
lation  of  the  £2  contribution  occurs  for  parabolic  bands. 

With  regard  to  the  third  moment  equation,  the  equation  for  energy  balance,  we 
find  nonparabollc  contributions  arising  even  for  uniform  fields.  First  It  Is 
noted  that  there  are  nonparabollc  contributions  to  the  energy  scattering  In¬ 
tegrals.  Of  more  significance.  It  Is  a  fact  that  the  time  derivative  of  the 
mean  thermal  energy  of  the  electrons,  3[3/2nkBTg]/ 3t,  becomes 
3 t3/2nkpTgg(Tg)]/3t.  The  effect  In  sample  calculations  of  the  contri¬ 
bution  g(Tg)  Is  to  Increase  the  time  to  relaxation.  Under  nonunlform  field 
conditions,  the  kinetic  energy  contribution  undergoes  similar  alterations. 

There  are  additional  changes,  each  of  which  can  be  traced  to  the  Inclusion  of 
the  f^  terra  In  the  expansion  of  the  distribution  function.  Again,  under  para¬ 
bolic  conditions,  these  additional  terms  undergo  specific  cancellation.  In 
summary,  the  presence  of  nonparabollc  band  structure  significantly  alters  the 
form  of  the  governing  equations  for  transient  transport  within  the  framework  of 
the  moments  of  the  Boltzmann  transport  equation.  The  question  next  Is:  How  do 
these  changes  alter  our  picture  of  velocity  overshoot  In  semiconductor  devices? 

The  typical  way  In  which  transient  transport  problems  are  attacked  under  uni¬ 
form  field  conditions  Is  as  follows:  1)  A  set  of  band  structure  parameters  are 
established  which  yield  steady  state  velocity  versus  field  curves  that 


are  compared  either  to  experiment.  If  available,  or  other  calculations.  2) 

Once  an  established  set  of  parameters  Is  obtained,  the  time  dependent  calcula¬ 
tions  are  performed.  The  results  of  these  calculations  whether  through  Monte 
Carlo  procedures,  moment  procedures,  or  Iterative  procedures,  all  show  similar 
velocity  overshoot.  The  wrinkle  In  this  procedure  Is  as  follows:  For  nonpara- 
bollc  bands  another  set  of  band  parameters  Is  chosen  from  which  a  mean  momentum 
versus  field  relationship  calculated.  From  an  equation  relating  velocity  to 
momentum,  the  resulting  steady  state  velocity  versus  field  Is  obtained.  Now, 
when  the  transient  calculations  are  performed,  the  transients  Involve  the  ap¬ 
proach  of  the  mean  momentum  to  steady  state,  from  which  the  velocity  transient 
Is  obtained.  From  a  practical  point  of  view.  It  will  not  surprise  anyone  that 
the  time  to  steady  state  with  and  without  nonparabollc  contributions  Is  dif¬ 
ferent.  But  we  are  not  Interested  solely  In  transients  to  steady  state.  For 
example:  When  a  negative  differential  mobility  element  Is  designed  to  operate 
as  an  active  device  under  94GHz  constraints.  It  Is  necessary  to  transfer  car¬ 
riers  from  the  T  to  L  and  X  valleys,  and  back.  One  criteria  for  sustained 
oscillations  Is  that  the  peak  to  valley  velocity  ratio  be  sufficiently  high. 
This  Implies  that  a  significant  amount  of  electron  transfer  from  the  T  valley 
has  occurred.  For  nonparabollc  bands.  Insofar  as  velocity  and  moment  are  not 
linearly  related,  there  are  situations  where  a  large  peak  to  valley  ratio  In 
velocity  does  not  necessarily  Imply  a  similarly  large  peak  to  valley  ratio  In 
momentum.  This  has  emerged  from  calculations  In  which  the  compound  semicon¬ 
ductors  were  placed  In  a  circuit  containing  reactive  elements.  Thus  a  new 
criteria  must  be  established  for  determining  the  upper  frequency  limit  for 
active  device  operation  In  which  nonparabollc  contributions  are  Introduced. 
Consider  next,  the  effects  associated  with  nonparabollclty  In  semiconductors. 

Calculations  were  performed  to  determine  specific  nonparabollc  effects.  These 
nonparabollc  effects  were  observed  for  GaAs,  InGaAs,  and  AlGaAs.  For  example: 
Steady  state  calculations  were  performed  for  GaAs  subject  to  parabolic  bands. 
The  field  dependence  of  velocity,  carrier  density,  and  electron  temperature  for 
each  of  the  three  valleys  were  obtained.  For  the  same  set  of  band  structure 
parameters,  e.g.,  phonon  frequencies.  Intervalley  deformation  potential  cou¬ 
pling  coefficients,  etc.,  the  calculation  was  redone,  but  with  a  nonparabollc 
dispersion  relationship.  Under  steady  state  conditions,  the  only  alteration  In 
the  equations  Is  In  the  scattering  Integrals.  The  results  of  the  calculation 
were  Initially,  partly  unexpected.  First,  as  anticipated  the  presence  of 
nonparabollc  bands  reduced  the  low  field  mobility  of  GaAs  from  Its  parabolic 
value.  The  unusual  result  was  that  the  alteration  In  the  scattering  rates  also 
effected  the  rate  of  electron  transfer.  The  result  being  that  more  carriers 
were  retained  In  the  F  valley  than  were  expected  from  parabolic  calculations. 
There  was  thus,  an  apparent  'Improvement'  In  the  'calculated'  velocity  field 
curve.  There  cannot,  of  course,  be  any  'Improvement'  In  the  velocity  field 
curve,  and  so  the  band  structure  parameters  for  the  GaAs  required  readjustment 
to  agree  with  experiment.  It  Is  noted  that  very  similar  results  were  reported 
by  A.  Sher,  at  the  recent  1986  DARPA  EHF  Review.  These  characteristics  were 
also  seen  for  ternary  materials. 

The  above  discussion  with  concerned  specifically  with  the  effects  of  nonpara¬ 
bollc  behavior.  The  InGaAs  and  AlGaAs  studies  were  performed  within  the 
context  of  broader  goals.  First,  It  was  anticipated  that  the  band  structure  of 
InGaAs  was  such  that  a  larger  fraction  of  carriers  could  be  retained  In  the 
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r  valley  than  was  possible  for  GaAs,  and  that  with  Its  high  mobility.  It  would 
offer  superior  current  drive  capabilities  for  submlcron  applications  than  GaAs, 
but  that  Its  Inferior  high  field  saturated  drift  velocity  would  make  It  less 
useful  than  GaAs.  Nonuniform  field  calculations  for  submlcron  GaAs  and  InGaAs 
were  performed  and  demonstrate  the  superior  current  drive  characteristics  of 
the  InGaAs.  An  additional  question  was  asked  with  respect  to  GaAs  versus 
InGaAs  (note.  In  this  discussion  we  are  concerned  with  a  50t  In  concentration); 
namely,  what  are  the  upper  frequency  limits  of  GaAs  versus  InGaAs.  Others  have 
addressed  this  question,  and  have  concluded  that  the  upper  frequency  limit  of 
GaAs  Is  higher  than  that  of  InGaAs.  Indeed,  a  recently  completed  DARPA  study 
performed  by  the  Principal  Investigator,  using  very  general  scaling  concepts 
corroborates  this  result.  The  final  conclusions,  however,  are  not  In.  For 
example:  Nonuniform  field  calculations  were  performed  through  solutions  to  the 
Boltzmann  transport  equation  to  simulate  Gunn  oscillations.  These  oscillations 
Involve  the  movement  of  dipole  layers  from  a  nucleatlon  site  to  a  collecting 
site.  Oscillations  In  GaAs  have  been  measured  at  frequencies  somewhat  In 
excess  of  lOOGHz.  Comparative  GaAs  and  InGaAs  calculations  were  performed. 

Both  of  them  require  different  levels  of  background  doping  to  achieve  high 
frequency  oscillations.  Both  semiconductors  displayed  oscillations  In  the 
lOOGHz  regime.  A  comparative  upper  frequency  limit  has  not  been  obtained  yet. 
Further  work  Is  needed.  It  Is  noted  that  In  the  nonuniform  field  calculations 
all  transport  was  assumed  to  be  parabolic. 

The  situation  with  AlGaAs  was  also  considered.  Here  calculations  were  perform¬ 
ed  with  varying  compositions  of  AlGaAs.  The  experimental  observation  of  a 
significant  dip  In  the  mobility  near  a  50%  composition  of  A1  was  not  observed, 
even  when  alloy  scattering  was  Included.  It  may  be  anticipated  that  such 
quantities  as  the  deformation  potential  coupling  coefficient  will  require 
significant  nonlinear  Increases  near  this  composition  range  to  account  for  the 
dip  In  mobility.  Major  attention  was,  however,  focused  on  a  30%  composition  of 
Al.  Within  this  range  AlGaAs  exhibited  a  marginal  region  of  negative  differen¬ 
tial  mobility.  The  Interest  In  the  AlGaAs  Is  clear  from  the  point  of  view  of 
the  new  heterostructure  studies,  and  calculations  were  performed  for  AlGaAs/ 
GaAs  HEMT  devices. 

The  HEMT  studies  were  performed  using  the  semiconductor  drift  and  diffusion 
equations  and  reveal  the  following  results,  some  of  which  were  Initially  sur¬ 
prising.  The  first  result  was  that  while  a  two  dimensional  electron  gas  formed 
at  the  AlGaAs/GaAs  Interface,  It  wasn't  as  had  been  expected  confined  to  a  nar¬ 
row,  sub-lOOA  width.  Rather,  there  was  significant  charge  Injection  within  the 
seml-insulatlng  substrate,  and  within  the  heavily  doped  AlGaAs  region.  This 
result  was  Initially  quite  controversial.  And  It  wasn't  until  very  recently, 
that  we  understood  the  origin  of  the  result  to  rest  with  the  dlffuslvlty  as¬ 
sumption  of  the  calculation.  For  example:  In  virtually  all  HEMT  studies,  both 
analytical  and  numerical,  the  Einstein  relationship  has  been  assumed.  This 
relationship  underestimates  the  high  field  dlffuslvlty  and  as  our  studies  In¬ 
dicated,  result  In  narrow  confinement.  We  have  used  an  empirical  relationship 
that  more  closely  represents  experimental  observation  and  permits  higher 
dlffuslvlty.  This  latter  result  Is  responsible  for  the  numerical  results. 

The  above  discussion  Involves  solutions  to  both  the  semiconductor  drift  and 
diffusion  equations  and  the  moments  of  the  Boltzmann  transport  equation.  As 
Is  well  known  there  are  difficulties  when  one  Is  considering  transport  at 


heterostructure  Interfaces.  To  overcome  these  difficulties,  we  have  been 
directly  Involved  In  formulating  quantum  transport  for  submicron  devices. 
Initially,  In  1981,  through  application  of  the  Wigner  distribution  function, 
a  series  of  moment  equations  were  formulated  for  pure  state.  In  the  absence  of 
scattering,  only  the  first  two  of  those  equations  were  independent.  These  two 
equations  were  equivalent  to  Schrodlnger’s  equation.  During  the  course  of  the 
ARO  study,  the  moment  equations  were  generalized  to  a  mixed  state.  In  general¬ 
izing  these  equations,  the  formulation  was  performed  In  terms  of  the  density 
matrix  rather  than  the  Wigner  function,  although  the  two  are  related  through  a 
transformation.  In  this  formulation  It  was  Immediately  apparent  that  a  separ¬ 
ation  could  be  established  Isolating  specific  quantum  contributions  from  the 
classical  contributions.  Further,  It  was  also  Immediately  apparent  that  no 
assumptions  on  the  properties  of  the  potential  were  necessary,  other  than  the 
fact  that  the  potential  seen  by  the  carriers  be  continuous. 

Each  of  the  above  studies:  Boltzmann  transport,  drift  and  diffusion,  and 
quantum  transport  have  been  presented  at  major  technical  meetings,  summer 
schools,  and  government  reviews.  They  are  all  being  prepared  for  technical 
publication.  Copies  of  publications  arising  from  this  study  are  Included  with 
this  report.  A  list  of  open  presentations  Is  given  below. 
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The  role  of  boundaries  and  interfaces  on  the  electrical  characteristics  of  lon^  and  siibnucuiu 
scale  compound  semiconductor  devices  is  discussed. 


1.  Introduction 

Theoretical  and  experimental  studies  of  compound  semiconductor  dovice.s 
over  the  past  decade  have  demonstrated  that  conditions  at  the  device  boundaries 
are  the  most  important  determinant  of  the  operating  characteristics  of  tlie 
device.  Studies  in  long  two-terminal  gallium  arsenide  devices  [1].  indium 
phosphide  (2],  germanium  [3],  and  cadmium  sulphide  (4]  have  demonstrated 
that  conditions  at  or  near  the  contacts  control  the  current-voltage  relations 
and  the  electrical  characteristics  of  any  resulting  instabilities.  For  two-terminal 
devices  it  is  often  possible  to  correlate  the  pre-  and  post-threshold  device 
behavior.  For  three-terminal  devices,  most  interface  studies  have  focused  on 
the  role  of  the  layer  just  under  the  principle  region  of  electron  transport.  In  the 
case  of  gallium  arsenide  field  effect  transistors,  heterostructure  interfaces  have 
been  incorporated  with  the  object,  e.g.,  of  confining  carriers  to  the  active 
region  in  the  case  of  a  heavily  doped  active  region,  [5,6],  or  of  providing  a  sea 
of  carriers  in  a  nominally  undoped  region. 

In  near  and  submicron  scale  devices,  the  relative  importance  of  the 
boundaries  increases  as  these  regions  occupy  a  si/.eable  fraction  of  the  device 
active  region  length,  and  the  up-  and  downstream  boundaries  begin  to  com¬ 
municate  with  each  other.  Controlling  these  boundaries  is  likely  to  be  the  most 
significant  task  of  device  physics  studies  in  the  immediate  future.  The  motiva¬ 
tion  for  this  outlook  is  based  upon  the  fact  that  transport  on  a  near  and 
submicron  scale  involves  nonequilibrium  effects  on  a  picosecond  lime  frame, 
and  its  effectiveness  is  based  on  transport  by  high  mobility  carriers.  I  hus, 
questions  such  as:  “How  many  carriers  are  injected  into  the  /'  valley  of  gallium 
arsenide,  and  with  what  energv  and  velocity?"  enter  the  picture  prominantly. 

()()VM.()2H/K.V000()  0()0()/.S().?.00  '  19S.?  Norih-Holland 


■  »•-  . 


V.*- 


H.L  iiruhtn.  J  P.  Krvskoxwkv  /  Ht^h  xpved  compound  senuiumtltu  for  tlcvu  vs 

The  purpose  i>f  tliis  paper  is  to  higlilight  the  effects  of  boundaries  on 
compound  semiconductor  devices,  particularly  gallium  arsenide.  In  doing  so, 
the  discussion  is  separated  into  two  parts;  (1)  the  role  of  boundaries  on  long, 
low  frequency  (<  20  (JHz)  gallium  arsenide  devices,  and  (2)  the  role  of  the 
boundary  on  high  frequency  near  and  submicron  scale  devices. 

For  long  devices  in  which  a  rich  body  of  experimental  work  exists,  a  review 
of  boundary  effects  is  given  in  terms  of  a  “pinned"  cathode  field  and 
“pinned"  cathode  current  model.  The  basis  for  these  models  is  the  notion  that 
contact  boundaries  may  be  phenomenologically  represented  as  either  tunneling 
or  thermionic  emi.ssion  dominated  regions,  with  a  varying  barrier  height.  The 
description  of  long  devices  assumes  that  all  pico.second  scale  transients  have 
occurred,  and  that  all  band  structure  population  statistics  are  adquately 
represented  by  steady  state  conditions. 

For  the  shorter  devices  it  is  the  short  time  transients  and  the  spatial  and 
temporal  transients  within  the  bands  that  are  calculated.  The  de.scription 
requires  solutions  of  the  Boltzmann  transport  equation  and  resolution  on  the 
scale  of  a  fraction  of  a  bulk  mean  free  path  is  needed.  These  problems  are 
discussed  in  section  3  where  a  review  of  such  phenomena  as  velocity  overshoot 
IS  given.  Fhe  ability  to  attain  contact  and  boundary  effects  permitting  the 
realization  of  the  high  overshoot  speeds  in  devices  is  the  core  on  which  the 
studies  of  .section  3  are  ba.sed. 

All  of  the  theoretical  studies  presented  arise  from  solutions  of  differential 
equations.  For  long  devices,the  “drift  and  diffusion"  equations  are  solved.  For 
short  devices,  moments  of  the  Boltzmann  transport  equation  are  solved.  In 
both  ca.ses  true  contacts  are  represented  as  boundary  conditions  to  the  dif¬ 
ferential  equations.  Thus,  the  studies  illustrate  the  effects  of  the  contacts 
and/or  interfaces  on  device  behavior.  The  results  of  these  studies,  when 
successful,  tend  to  highlight  what  is  unknown  about  device  material  parame¬ 
ters.  For  the  case  of  long  devices,  it  is  the  cathode  boundary  condition  that  is 
unknown.  Here  the  sensitivity  of  the  results  to  varying  the  cathode  field 
identifies  the  significance  of  the  contacts  on  device  behavior  1 1  ].  Similar 
sensitivity  lUiidies  are  discussed  in  connection  with  solutions  of  the  Boltzmann 
transport  equation,  where  a  range  of  parameters  is  chosen  to  identify  the 
conditions  for  "injecting"  and  “blocking"  contacts.  A  sensitivity  analysis  of 
the  electronic  contribution  to  the  thermal  conductivity  is  also  included  as  its 
effects  are  dramatic. 


2.  Boundary  conditions  to  negative  differential  conductivity  devices 

In  this  section  a  brief  review  of  the  influence  of  boundaries  on  the  behavior 
of  NDC  devices  is  given.  Fig.  I  displays  typical  boundary-dependent  data  from 
three  different  gallium  arsenide  two-terminal  devices  (7).  'Fhe  lower  portion  of 
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each  diagram  displays  current  versus  voltage  characteristics,  wliile  the  upper 
portion  shows  voltage  versus  distance  at  one  bias  point.  I'ig.  la  shows 
metisurements  for  a  device  in  which  the  metal  c<intacts  are  far  removed  frtmi 
the  active  region  of  the  device.  The  current- voltage  rehition  is  relatively  linear 
until  a  point  where  current  oscillations  occur  I  he  field  profile  just  prior  to  the 
oscillation  is  relatively  uniform  within  the  active  region  of  the  device,  and  is 
near  zero  at  the  ends  of  the  active  region.  Fig.  lb  represents  a  set  erf 
metisurements  in  which  the  metal  contact  abutted  the  active  regiern  of  the 
device.  I  he  current- voltage  characteristics  remained  linear  to  threshold  which 
again  was  manifested  by  a  current  oscillation.  Notably  different  here  is  the 
lower  average  field  prior  t(r  the  instability  and  the  enhanced  voltage  drop  at  the 
cathode,  big.  Ic  displays  results  for  another  device  with  a  metal  contact 
abutting  the  active  region,  lor  this  case  there  is  a  sublinear  current  virltage 
characteristic  and  no  instability.  I  he  probed  yoltage  versus  liistance  shows  a 
large  voltage  drop  at  the  vicinity  of  the  cathode. 

The  electrical  characteristics  associated  w'ith  fig.  I  have  been  describcti  as 
representative  of  “ohmic”  (fig.  la),  “slightly  blocking”  (fig.  lb),  and  “strongly 
blocking"  (fig.  Ic)  contacts.  One  of  the  earliest  models  employed  for  c.yplaimng 
the.se  results  assumed  a  “pinned”  value  of  cathode  electric  field  (1).  Other 
models  in  which  the  cathode  conductivity  [8]  or  doping  [9]  have  also  been 
suggested  w'ith  varying  degrees  of  success.  The  “pinned”  cathode  field  model 
developed,  partially  as  a  consequence  of  the  w-ay  the  governing  equations 
de.scribing  current  instabilities  was  written.  Here  the  one-dimensional  differen¬ 
tial  equation  for  total  current, 

,  .  /  ihv  \  a/: 

was  rewritten,  using  Poisson's  equation,  as 


a-/-;  a/: 


This  is  a  second  order  partial  differential  equation  requiring  t\so  boundary 
conditions  on  /.'(  .V,  7  )  and  one  initial  condition.  In  the  above  7(  7  )  represents 
current  density.  A'  is  carrier  density,  F  and  I)  are  field  dependent  selocity  and 
diffusion. 

Srrlutions  ttr  eq.  (2)  have  been  used  successfully  to  simulate  device  results 
similar  to  those  of  fig.  1.  F'or  gallium  arsenide  in  which  the  threshold  field  for 
negative  differential  nuibilitv  is  appro.ximately  .a. 2  kV/cm  =  /-  in,  qualitatively 
similar  results  occur  for  solutions  with  cathode  fields  falling  into  any  one  of 
the  following  three  groups:  0  s,  /:(  X  =  0,  7  )  /: , ,,  ^  /:'(  .V  =  0,  /')  •-  4 

/:  I /•-'(  ,\  ^  0,  I  )  >  4  The  simulations  with  pinned  fields  falling  m  either 
group  I.  2  or  }  yiekl  electrical  characteristics  similar  to  those  of  figs  la,  lb  ami 
Ic.  respectively.  The  crucial  fe.iture  of  this  model  is  th.it  the  cathode  field  is 
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and  cathode  accumulation.  From  ref.  (10],  with  permission. 
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pinned,  necessitating  that  any  instabilities  in  current  occur  at  a  critical  value  of 
current  density.  The  field  profiles  a.ss<Kiated  with  cathode  fields  in  the  range 
0  <  £■(  ,V  =  0,  T)  <  and  /: , ,,  <  /:(  A'  =  0.  T)  <  4  /: , ,,  are  sketched  in  figs. 
2a  and  2b  re.spectively.  For  reference,  a  velocity  field  curve  with  velocity  scaled 
to  current  as  qN^y(E)  =  and  with  a  region  of  negative  differential 

mobility  is  also  included.  Fig.  2  is  understood  as  follows:  The  second  column 
of  each  section  shows  the  electric  field  versus  distance  profiles.  E(  X)  begins 
with  a  value  E^  at  the  cathode  and  extends  downstream  to  a  value  By 

current  continuity,  the  current  everywhere  within  the  device  is  given  by 
J  =  qNV(Ey^).  For  J  <  E^),  a  region  of  charge  depletion  forms  near  the 

cathode  for  E^  below  and  within  the  region  of  negative  differential  mobility 
(NDM).  Increasing  the  current  until  J  =  J„(E^)  introduces  charge  neutrality 
everywhere  for  E^  <  Ey^^.  However,  because  E  is  a  double  valued  function  of 
V,  for  E^  within  the  NDM  region,  approximate  charge  neutrality  exists  near 
the  cathode  for  J  <  £  ),  and  for  regions  .sufficiently  far  downstream  from 

the  cathode.  Charge  neutrality  breaks  down  between  these  two  regions.  Finally. 
{or  J  >  J^(E^)  an  accumulation  layer  forms  near  the  cathode.  For  fig.  2a.  the 
accumulation  layer  is  stable  until  the  bulk  field  exceeds  £th-  ^b.  the 

accumulation  layer,  followed  downstream  by  the  depletion  layer,  is  often 
unstable  and  leads  to  cathode  originated  instabilities  (10). 

The  situation  corresponding  to  fig.  Ic  is  often  represented  by  very  high 
cathode  fields.  The  field  profiles  are  tho.se  appropriate  to  a  wide  region  of 
charge  depletion  near  the  cathode.  The  profiles  are  electrically  stable. 

The  characteristic  feature  of  these  nonuniform  field  profiles  is  that  their 
structure  is  significantly  affected  by  the  field  being  a  double  valued  function  of 
velocity.  The  pinning  of  the  cathode  field  is  not  necessarily  common,  however, 
to  all  semiconductor  devices.  For  example,  it  was  also  applied  to  InP  devices, 
where  it  w'orked  for  a  significant  number  of  cases.  However,  a  broad  class  of 
device  behavior  could  not  be  accounted  for  through  its  u.se  [2).  These  devices 
showed  anomalously  high  efficiency  and  significantly  low  DC  current  Icsels. 
Spontaneous  Gunn  type  oscillations  did  not  occur.  Rather,  device  operation 
required  a  tuned  circuit.  The  details  of  the  o.scillalion  were  thought  to  depend 
critically  on  the  cathode  boundary  condition,  which  in  this  case  was  taken  as  a 
fixed  cathode  conduction  current  (10,1 1]. 

The  distinction  between  “pinned”  cathode  field  and  “pinned"  cathode 
conduction  current  is  placed  in  perspective  in  fig.  3  and  in  the  following 
equation 

y(7')=y..(£  )+e^.  (3) 

Flq.  (3)  is  the  equation  for  total  current  through  the  boundary  to  the  device, 
y  (£\  )  represents  the  current-field  relation  at  the  cathode  (121  which  mav  be 
expected  to  differ  from  that  of  the  semiconductor  device.  Two  such  tvpcs  of 
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l  ig.  3.  Cathode  conduction  current  curve  for  an  approximately  "pinned"  cathode  field,  curve  A; 
and  pinned  cathode  conduction  current  (at  high  field  fields),  curve  B.  Curves  are  obtained  from  eq. 
(4).  Also  shown,  for  reference,  is  the  neutral  current  versus  field  curve  for  GaAs.  From  ref.  (lO), 
with  permission. 


curves  are  represented  by  curve  A  and  curve  B  of  fig.  3.  Curve  A  is  closely 
related  to  the  pinned  cathode  field  model  while  curve  B  is  associated  with  the 
pinned  cathode  current  model.  The  similarity  in  “form”  of  curves  A  and  B  to 
moderate  barrier  height  tunneling  and  thermionic  emission  dominated  con¬ 
tacts  is  deliberate  (10),  and  the  equation  used  to  arrive  at  these  curves  is  shown 
below 


J.  (  C  )  =  -  V, 


nkr 


-  -  1 
n 


k  T 


which  was  adapted  from  studies  on  the  unalloyed  mctal/semiconductor  con¬ 
tact  (13).  Its  use  here  presumes  a  similar  description.  For  the  unalloyed  contact, 
n  IS  the  ideality  factor  and  describes  the  contact  as  dominated  by  thermionic 
emission  (n  =  1)  or  by  tunneling  (n  »  I).  is  the  reverse  current  flux  and 
may  be  related  to  the  barrier  height  phenomenologically  through  the  Richard¬ 
son  equation  (10). 

Detecting  a  particular  contact  effect  on  a  device  is  a  difficult  procedure.  For 
long  devices  current  voltage  characteristics  as  represented  by  fig.  1  are  often 
signatures  of  a  contact  classification.  For  short  devices  proximity  effects 
introduce  an  additional  complication  and  current-voltage  measurements  are 
less  valuable.  One  type  of  measurement  which  may  serve  to  provide  informa¬ 
tion  about  the  boundary  is  a  mrise  measurement. 

Here  the  situation  to  envision  is  that  if  the  field  is  pinned  within  the 
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negative  differential  mobility  region,  increasing  the  bias  will  result  in  an 
increase  in  the  length  of  the  NDM  region.  Any  fluctuation  will  sustain 
enhanced  amplification  and  the  noise  will  increase.  If  increasing  the  bias 
resulted  in  carrier  injection  into  the  device,  the  field  at  the  cathode  is  likely  to 
decrease  with  increasing  bias  and  the  noise  is  expected  to  decrea.se.  While  these 
results  should  be  folded  in  with  the  field  dependence  of  velocity  and  diffusion, 
a  simple  analytical  noi.se  calculation  a.ssuming  a  three  piece  linear  approxima¬ 
tion  to  represent  GaAs  has  been  performed  |14j. 

In  this  calculation,  the  “impedance  field  method"  1 15]  is  applied  to  calculat¬ 
ing  noise  due  to  thermal  velocity  fluctuations  amplified  within  the  device.  The 
mean  .squared  noise  voltage  per  unit  band  width  (15)  is  computed. 

=  |VZpN/ni(vol).  (5) 

where  vZ  is  the  impedance  field  vector  ( I5j.  I  he  calculation  is  performed  for  a 
10  /xm  long  element  with  a  doping  of  l()'‘’/cm’.  The  element  sustains  the  field 
profile  of  fig.  2b  where  it  is  seen  that  the  NDM  region  increases  with 
increasing  bias.  The  calculations,  which  are  discus.sed  in  detail  eLsewhere.  [14] 
are  expressed  in  terms  of  the  noise  figure  ( I6|: 


NF=  1  -I- 


1 

4/cot1«|' 


(6) 


where  R  is  the  real  part  of  the  device  impedance.  The  results  of  the  calculation 
are  displayed  in  fig.  4,  where  the  noi.se  figure  is  sketched  as  a  function  of  bias 
current  and  transit  angle.  The  results  appear  as  a  signature  of  the  effects  of  the 
cathode  boundary.  First,  at  low  values  of  transit  angle  9  =  u7'{A).  where  T(A) 
is  the  transit  time  across  the  negative  differential  mobility  region,  the  noise 
figure  increases  with  increasing  bias.  Ibis  corresponds  to  an  increase  in  the 
length  of  the  negative  differential  mobility  region  and  enhances  amplification 
of  any  fluctuation  originating  there.  More  interesting  structure  is  present  at 
higher  frequencies  and  higher  bias  where  the  noise  figure  increases  and  then 
shows  a  singularity.  On  the  other  side  of  the  singularity  there  is  a  “U’'-shaped 
region  ending  again  at  a  singularitv.  The  strong  increa.se  in  noi.se  figure 
repre.sents  the  approach  of  |/?|  =>  0.  Here,  at  low  frequencies,  the  real  part  of 
the  impedance  is  positive,  and  becomes  negative  at  frequencies  somewhere 
between  -n  ^9  ^  2tt.  In  going  from  positive  to  negative  values  it  passes  through 
zero,  hence  the  singularity.  The  frequency  range  for  small  signal  negative 
resistance  increa.ses  with  increasing  bias  |I0]  reflecting  the  broadening  of  the 
negative  differential  mobility  regions  -  a  broad  “U"-shaped  region  appears. 
Both  the  increasing  noi.se  figure  at  low  frequencies,  and  the  “U"-.shaped  region 
at  high  frequencies  are  characteristics  of  an  increasing  depletion  layer  uidth. 
Note  that  increasing  the  bias  still  further  will  result  in  an  electrical  instabilitv 
(10). 
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transit  angle,  B 

Fig.  4.  Noise  figure  versus  transit  angle  and  normalized  bias  current  j  =  -/h/A'ocFj,,  for  a  10  fim 
long  element  with  a  depletion  layer  profile.  From  ref.  (14),  with  permission. 


3.  Klectron  transport  in  near  and  submicron  devices:  the  role  of  the  boundaries 

The  discussion  of  the  above  sections  dealt  with  devices  whose  lengths  were 
typically  10  /am  long  or  longer.  Transport  for  these  devices  is  generally 
discussed  though  use  of  the  drift  and  diffusion  equation  (eq.  ( 1)).  For  near  and 
submicron  length  devices  electron  transfer  is  generally  not  complete  until  a 
substantial  fraction  of  the  device  has  been  traversed.  Consequently  the  drift 
and  diffusion  description  is  inadequate  and  solutions  to  the  Boltzmann  trans¬ 
port  equation  are  required. 

F-or  two  level  transfer  in  GaAs  the  steady  state  velocity  field  curve  for 
carriers  in  the  centra!  and  satellite  valleys  are  shown  in  fig.  5,  where  the  net 
velocity  V  is 

F= /V,F,  +  yV.K,.  (7) 


Here  the  subscript  1  denotes  transport  in  the  central  valley  (/'  valley  for  GaAs) 
and  the  subscript  2  denotes  transport  in  the  satellite  valley.  Most  device  design 
is  concerned  with  controlling  the  time  spent  by  the  numbers  of  carriers  m 
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in  valleys  1  or  2;  and  circuits,  interfaces  or  contacts  are  sought  which  will  allow 
for  this  control.  Solutions  to  the  Boltzmann  transport  equation  which  provide 
the  required  nonequilihrium  transport  behavior  are  obtained  by  a  number  of 
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big.  5.  Steady  state  velocity  versus  field  curves  for  /'  and  satellite  valleys  for  gallium  arsenide 
From  ref.  (22],  with  permission. 


different  techniques,  the  most  familiar  t)f  which  is  the  Monte  Carlo  method 
[17].  Other  methods  include  the  Rees  iterative  technique  [18]  and  balance 
equation  solutions  [19],  which  are  considered  below. 

The  balance  equations  discussed  below  are  the  first  three  moments  of  the 
Boltzmann  transport  equation.  They  form  a  subset  of  the  infinite  hierarchy  of 
moment  equations,  and  as  such  do  not  form  a  clo.sed  set  of  equations.  F.ach 
equation  introduces  a  higher  order  moment  not  defined  by  the  given  set  of 
equations.  The  most  common  form  of  the  balance  equations  is  obtained  by 
assuming  a  distribution  function  of  the  form 


/  =  /),  “-TT  exp 
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for  each  species  More  generally  [20] 

/  a  a’  a’  \ 

.w;  +  '.wju,  *  7w. 


where  it  is  insisted  that  the  first  term  in  the  expansion  gives  the  correct  local 
values  of  density,  velocity  and  energy.  The  coefficients  in  the  above  expansion 
are  inversely  proportional  to  density  and  are  model  dependent.  The  first 
correction  to  the  local  equilibrium  balance  equtions,  the  so-called  hydrody¬ 
namic  approximation  yields 
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wtiere  t  is  an  electron  temperature. 
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a  is  a  stress  tensor  arising  from  nonuniform  velocilv  distributions,  and  k  is  the 
thermal  conductivity.  The  stress  term  is  dissipative  in  that  when  a  nonuniform 
velocity  distribution  is  impressed  on  an  electron  stream,  there  will  be  reactive 
forces  tending  to  smooth  them  out.  At  this  pt>int.  lunvever.  these  terms  are 
regarded  as  pheiumienological  entries. 

d'he  boxesl  terms  ab<ive  represent  ct>ntributions  from  the  nonspherical 
nature  of  the  distribution  function.  I  hese  contributions  liase  not  been  included 
in  evaluating  the  collision  integrals,  i.e  ,  (  .  )  ^  hut  uill  be  discussed  in  a 
future  paper.  I  he  coIIiskui  integrals  are  discusscil  m  lel  |2lj  I  he  underlined 
and  boxed  terms  are  ignored  m  the  ilrifi  .iiul  liiffusu'ii  .ippioximation. 

Under  uniform  held  conditions  the  ine.in  response  i>|  i.imers  to  a  sudden 
change  m  electric  field  is  shown  m  fie  ('  I  he  high  pe.ik  velocits  is  a 
consequence  of  carriers  being  ret.imed  m  the  I  \.ille\s  for  time  rlui.ition 
upwards  ot  0.2  ps  before  uiulergoing  ir.msfei  to  the  upper  v.illevs  The  peak 
velrK'itv  111  fig.  (i  IS  extremelv  high  It  lepiesenis  .ui  upt>ei  bouiui  on  the  c.irner 
velocitv  th.it  m.i\  be  expecieii  .n  this  lield  .uul  pi. Aides  the  nu'tis.ition  for 
designing  bound. ines  tli.it  peimil  .uliiesemeni  o|  these  \  .dues 

I  he  pe.ik  selocits.  hovsesei.  is  seiisiuse  lo  t>i.is  nse  lime  .md  hints  at 
problems  to  be  .iddiessed  in  ilesigiiiiu'  ippi  opi  i.iie  divue  houiul.iis  condi¬ 
tions  .See  (ig  ^  I  tnlheiinofe  bia.iiise  i  Innli-  liiiii-  i.  leqiiiied  loi  c. liners  to 
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Fig.  6.  Transient  velocity  versus  time,  for  electrons  subject  to  a  sudden  chinge  in  field  of 
magnitude  27  kV/cin.  From  ref.  [22|.  with  permission. 


transfer  into  subsidiary  portions  of  the  conduction  band,  on  a  near  and 
submicron  scale,  the  velocity  field  characteristics  arc  expected  to  be  length 
dependent.  See  fig.  8.  Thus,  near  and  submicron  devices  will  categorically  be 
sensitive  to  both  boundary  effects  and  device  length,  neither  of  which  will  be 
independent  of  the  other. 

Figs.  5  to  8  are  signatures  of  high  speed  submicron  transport.  What  are  the 
consequences  for  device  behavior?  Can  a  suitable  .set  of  contact  or  interface 
conditions  be  found  to  achieve  high  speeds  associated  with  overshoot?  While  a 
detailed  study  on  a  scale  similar  to  that  for  long  devices  is  not  yet  available, 
some  results  are  known.  We  discuss  these  below  for  a  collection  of  2  /xm  long 
devices,  essentially  undoped.  A',,  =  5  X  lO'^/cnv’.  each  subjected  to  a  bias  of  2 
V.  Thc.se  one-dimensional  studies  are  in  steady  state,  i.e..  d/3T=0.  and  are 
subjected  to  the  following  cathode  boundary  conditions: 


60K 
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For  A  positive  (iicgativo)  local  charge  accumulation  (depletion)  occurs  at  the 
cathode  boundary.  B  is  generally  greater  than,  or  equal  to.  300  K  and  is  a 
measure  of  the  mean  thermal  energy  of  the  /'-valley  carriers.  As  in  section  2, 
the  results  are  placed  into  three  categories,  “ohmic”,  “slightly”  and  “strongly" 
blocking  contacts.  The  “ohmic”  results  are  shown  in  fig.  9. 

For  “ohmic”  boundaries  an  appropriate  set  of  constants,  with  ref.  1 19].  are: 
A  -  0.2.  =  12.000  enr/V  •  s.  and  t,  =  300.  For  this  case  carriers  enter  the 

device  with  speeds  greater  than  that  associated  with  the  central  valley  mobility. 
Consequently  there  is  an  accumulation  of  carriers  at  the  boundaries,  resulting 
in  low  values  of  cathode  field.  The  field  starts  off  at  nearly  1  kV/cm  and 
approximately  2500  A  must  be  traversed  before  significant  transfer  occurs. 
Increased  transfer  results  in  a  lowering  of  the  mean  carrier  velocity,  necessitat¬ 
ing  an  increase  in  mobile  charge  as  the  anode  is  approached.  The  average 
velocity  across  this  device  obtained  from  the  relation 

K.  =-^/‘7A'o  (17) 

is  =  1.78  X  10"  cm/s.  with  A,)  =  5  X  lO'^/cm^.  The  advantage  of  overshoot 
is  not  fulfilled  by  this  contact  device  length  configuration. 

Lowering  the  cathode  mobility  to  a  value  below  that  associated  with  the  /' 
valley  results  in  a  more  rapid  dispersal  of  carriers,  and  cathode  adjacent  charge 
depletion,  associated  with  slightly  blocking  contacts,  occurs.  This  is  seen  in  fig. 
10  for  =  -0.1 1,  =  6000  enr/V  •  ,s,  and  t,  =  300  K.  It  is  noticed  that  the 

cathode  field  for  this  case  is  approximately  4  kV /cm,  which  is  higher  than  that 
associated  with  the  “ohmic”  contact  condition  of  fig.  There  are,  however, 
important  similarities  between  figs.  9  and  10.  In  both  cases  the  carriers 
adjacent  to  the  cathode  are,  for  all  practical  cases,  F-valley  electrons.  Very  little 
transfer,  which  is  determined  by  carrier  energy  (temperature),  has  occurred.  In 
addition,  sufficiently  downstream  from  the  cathode  the  carriers  appear  to  be 
ignoring  the  cathode  condition  and  are  dominated  by  the  downstream  \oliage 
drop,  whose  spatial  distribution  is  about  the  same  for  both.  Tlie  average 
velocity  for  this  case  is  1',.  =  1.75  X  lO’  cni/s,  slightly  below  that  of  fig.  9. 

A  significant  change  occurs  when  the  mean  energy  of  the  I  -valley  entering 
carriers  is  elevated.  F-or  the  parameters  of  fig.  10.  but  with  t,  =  1200  K..  a 
substantial  amount  of  transfer  occurs  at  the  cathode,  resulting  in  a  lowering  of 
the  current  through  the  device.  The  cathode  field  is  approximately  7  kV/cm, 
higher  than  that  as.sociated  with  fig.  10.  and  the  downstream  field  is  lower  (see 
fig.  1 1).  The  average  velocity  of  the  carriers  in  this  case  is  lowered  to  1.28  x  10^ 
cm/s,  even  though  the  central  valley  carriers  are  traveling  at  higher  speeds  (sec- 
fig.  5). 

The  pre.sence  of  moderately  high  cathode  fields  is  attractive  if  a  sufficient 
number  of  carriers  can  be  retained  in  the  central  valley  where  they  can  sustain 
high  transit  vel(x:ities.  While  this  case  is  discussed  in  more  detail  below,  the 
simple  ruse  of  injecting  excess  carriers  into  a  device  with  the  contact  conditums 


ELECTRIC  FIELD  (kV/cm)  FRACTIONAL  POPULATION 
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Fig.  I  I.  As  in  fig.  9,  but  with  =  6000  cm^V  s.  A  =  -0.1  1  and  li  =  1200  K. 


when  compared  to  fig.  1 1.  In  this  case  =  1.30  x  10’  cm/s. 

The  situation  of  a  strongly  blocking  contact  is  illustrated  in  fig.  13  for 
A  =  9.0,  n  =  1000  crn^/V  ■  s,  and  T  =  3000  K.  For  this  ca.se  virtually  all  carriers 
are  swept  away  from  the  cathode,  with  the  satellite  valley  carriers  accounting 
for  most  of  the  transport  in  the  transition  region.  The  cathode  field  is 
extremely  high,  approaching  60  kV/cm.  The  average  velocity  is  approximately 
zero. 

The  above  results  focus  attention  on  the  role  of  the  upstream  boundary  in 
the  distribution  of  charge  near  the  interface.  The  results  are  more  general  than 
those  obtained  using  the  steady  state  field  dependent  velocity  relation,  espe¬ 
cially  in  identifying  the  fact  that  electron  transfer,  even  in  the  case  of  partially 
blocking  contacts,  may  not  occur  until  some  point  downstream  from  the 


('12  //./.  trrufiin.  J  /’  \k.»t  \k  \  Htsih  i  ttniptninii  \cinn  Dtiiiiu  tor  r/<’f 


Fig.  13.  As  in  fig.  9,  but  with  =  1000  cmV^  s.  -9.0  and  B  -  3000  K. 


cathode  contact.  The  implication  of  this  is  that  the  detailed  role  of  the  contact 
in  controlling  current  instabilities  may  need  generalization.  On  another  matter, 
it  is  apparent  that  the  advantages  of  velocity  overshoot  are  ab.sent  from  the  2 
nm  long  calculation.  The  reason  for  this  is  direct:  the  carriers  have  traversed  a 
path  of  sufficient  length  to  exceed  the  threshold  energy  for  electron  transfer.  In 
this  case  electron  transfer  results  in  velocity  saturation  and  currents  ap¬ 
propriate  to  the  steady  state  parameters.  Thus,  while  a  blocking  contact  may 
be  expected  to  yield  high  entrance  vekK'ities  for  one  species  of  carriers,  electron 
transfer  prevents  a  large  fraction  of  these  carriers  from  enjoying  the  speed 
advantage.  The  message  here  is  that  a  forgiving  contact  is  useful  only  when 
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coupled  to  an  appropriate  device  design.  Witliin  the  framework  of  the  dis¬ 
placed  Maxwellian,  improvements  in  device  structure  are  synonymous  with 
values  of  electron  temperature  that  are  below  that  required  for  electron 
transfer. 


Nq  =  5  X  lO'*  /cm* 


DISTANCE  (microns) 


3 


DISTANCE  (microns) 

Kig.  14.  (a)  Doping  distribution  for  N-N  -N  device.  Bias  is  I  V,  and  conditions  at  the  cathode 
boundary  are  as  in  fig.  (b)  F.lectric  field  distribution. 
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One  way  of  reducing  the  electron  temperature  is  simply  to  go  to  shorter 
device  structures,  and  recent  results  for  gallium  arsenide  at  room  temperature 
suggest  that  device  structures  not  greater  than  0.25  ;tm  may  be  required.  A 
structure  that  has  recently  been  discussed  within  the  framework  of  near  and 
submicron  devices  is  the  N-N  -N  structure  shown  in  fig.  14.  This  structure 
has  four  interfaces  to  contend  with.  The  most  significant  a.spect  of  this 
structure  is  that  it  introduces  an  abrupt  change  in  field  at  the  N-N  interface, 
which  tends  to  emphasize  overshoot  contributions.  The  detailed  input  to  this 
calculation  is  as  follows.  The  N  region  is  doped  to  5  X  lO'^/cm^  with  the  N" 
region  an  order  of  magnitude  less  in  doping.  For  this  device  the  cathode  ohmic 
conditions  of  fig.  9  were  imposed,  and  an  average  field  of  5  kV/cm  was 
imposed.  We  note  that  over  a  distance  of  approximately  1.0  ^m  there  is  very 
little  electron  transfer.  Current  is  carried  mainly  by  T-valley  electrons  who<:e 
carrier  velocity  peaks  near  6  x  10’  cm/s,  showing  substantial  overshoot.  The 
electron  temperature  for  this  calculation  is  shown  in  fig.  15,  and  is  reasonably 
low.  for  this  case  is  1.3  X  10’  cm/s. 

The  results  of  this  calculation  are  very  encouraging.  However,  they  are  not 
solely  a  consequence  of  structural  device  changes.  The  retention  of  /'-valley 
electrons  in  the  fig.  15  calculation  is  a  consequence  of  using  an  electronic 
thermal  conductivity  appropriate  to  5  X  through  the  Wiedemann- 

Franz  ratio.  A  reduction  in  electronic  thermal  conductivitiy  results  in  steeper 
gradients  in  electron  temperature.  Fig.  16  shows  results  for  a  value  of  thermal 
conductivity  in  which  substantial  transfer  occurs  within  0.5  fim  of  entry  into 
the  N'  region.  The  average  current  has  dropped  to  6.5  x  10*  cm/s.  The  first 
set  of  results  which  virtually  eliminate  transfer  must  be  regarded  as  an  upper 
bound  on  current,  whereas  the  second  set  is  a  lower  bound.  Actual  device 
results  are  likely  to  fall  between  the  two.  One  important  measurement  likely  to 
provide  significant  information  is  whether  a  1  jim,  10' Vcm’  element  will  show 
small  signal  gain. 

Returning  to  the  boundary  conditions,  there  are  several  points  to  be  made. 
First,  there  is  a  clear  indication  as  to  the  procedures  necessary  for  generating 
an  entire  set  of  “ohmic”  and  blocking  contacts.  It  is  not  likely,  however,  that 
the  boundary  prescriptions  of  figs.  9  to  13  are  unique;  other  sets  of  conditions 
can  be  envisioned  to  yield  similar  results  under  DC  conditions.  Distinguishing 
between  different  boundary  condition  effects  will  come  from  time  dependent 
studies.  The  question  then  is:  “How  model  dependent  are  the  results?” 

The  governing  equations  are  derived  starting  from  the  condition  of  a 
displaced  Maxwellian,  which  assumes  strong  electron-electron  interaction. 
Boundary  scattering  is  likely  to  substantially  alter  this  interaction  -  near  the 
boundary  [25j.  Results  that  appear  to  be  model  independent  are  tho.se  a.s.six;i- 
ated  with  the  entrance  velocity.  If  central  valley  carriers  enter  with  speeds 
greater  (less)  than  that  dictated  by  equilibrium  band  structure  considerations, 
carrier  accumulation  (depletion)  will  occur. 
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Fig.  17.  Represenlalion  of  sputial  dependence  of  current  flow  through  two  level  device  as  a  serial 
chain  of  parallel  elements.  The  first  subscript  in  each  element  represents  current  through  the 
central  valley  1  or  satellite  valley  2. 
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Another  point  worth  raising,  particularly  with  regard  to  the  absence  of 
overshoot  in  figs.  9  to  13,  is  the  fact  that  for  at  least  half  of  these  devices  the 
space  distributions  were  in  steady  state  equilibrium  with  results  similar  to  those 
of  the  drift  and  diffusion  equations.  While  going  to  shorter  device  lengths  does 
not  necessarily  prevent  this  from  occurring  (26],  the  length  dependence  has  not 
been  explored. 

The  results  of  figs.  9  to  13  also  suggest  that  transport  in  a  multivalleyed 
system  may  be  controlled  by  one  of  the  valleys.  In  the  above  examples,  control 
is  by  the  central  valley;  satellite  valley  boundary  conditions  are  relatively 
benign.  To  see  the  reasons  for  this  it  is  necessary  to  turn  to  fig.  17,  which 
represents  current  flow  through  the  nonlinear  elements  as  a  linear  chain  of 
parallel  resistors.  In  this  figure 

^  =  '=1.2 . N, 

4  =  A'2,^"2,.  '=1.2 . N. 

The  voltage  across  element  i  is  F,AX.  Thus  if,  e.g.,  the  F  valley  in  element  1 
sustains  a  net  carrier  density  below  the  background,  it  will  have  a  higher 
resistance  and  voltage  drop  than  an  element  with  a  density  closer  to  back¬ 
ground.  This,  of  course,  is  consistent  with  the  results  of  fig.  1 1  and  suggests 
that  the  portion  of  the  conduction  band  with  boundary  conditions  most 
strongly  departing  from  the  uniform  field  conditions  will  be  the  dominant 
boundary  condition. 

The  contact  boundary  conditions  discussed  in  figs.  9  to  13  address  only  a 
small  part  of  the  problem.  It  is  certainly  unrealistic  to  assume  that  undesirable 
departures  in  doping  will  be  absent.  A  simulation  of  a  nominally  10''/cm’ 
doped  device,  with  a  10%  decrease  in  doping  over  a  distance  of  1000  A,  is 
shown  in  fig.  18.  The  distortion  in  electric  field,  for  the  .set  of  boundary 
conditions  listed  in  the  caption,  is  such  as  to  prevent  any  real  overshoot  from 
occurring. 

The  above  results  which  reflect  the  influence  of  space  charge  on  transport  in 
devices  .should  be  compared  to  uniform  field  calculations  to  indicate  the  goals 
that  perhaps  should  be  taken  for  some  device  structures.  Fig.  19  displays 
velocity  versus  distance  curves  for  electrons  subjected  to  a  sudden  value  of 


fractional  population  field  (kV/cm) 


Fig.  20.  Velocity  versus  distance  for  electrons  with  varying  entrance  energies,  subject  to  a  sudden 
change  in  electric  field.  From  ref.  (27),  with  permission. 


4.  Summary 

The  experimental  situation  is  such  that,  with  the  exception  of  long  com¬ 
pound  semiconductor  devices,  there  are  very  few  data  on  the  role  of  boundaries 
and  contacts  to  submicron  devices.  The  reason  for  the  paucity  of  data  lies  in 
the  fact  that  most  submicron  devices  are  three  terminal  device  designs  and  the 
third  terminal  tends  to  mask  the  role  of  the  contact  boundaries.  This  is 
extremely  unfortunate  since  it  is  likely  that  two  terminal  device  measurements 
will  indicate  what  can  be  achieved  in  controlling  the  entrance  dynamics  of  the 
carriers.  To  date,  most  two  terminal  device  measurements  on  simple  device 
structures  have  concentrated  on  the  role  of  transport  within  the  device,  and 
rai.se  the  question  of  whether  “ballistic”  motion  is  possible  (28).  Ba.sed  on  the 
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liistory  of  vacuum  tuhc  ilyiiamics  (29|  it  should  be  recalled  that,  if  transport  is 
ballistic,  the  electrical  characteristics  will  be  controlled  by  the  contacts. 

riie  situtation  in  submicron  devices  is  further  complicated  by  communica¬ 
tion  between  the  up-  and  down-stream  contacts.  Thus  it  may  be  expected  that 
the  influence  of  a  bltK'king  contact  on  the  electrical  characteristics  of  long  and 
submicron  devices  will  be  different.  For  .submicron  devices  simple  current 
voltage  measurements  may  be  rendered  useless  as  a  diagnostic  tool.  This  is 
certainly  not  the  case  in  long  devices. 

The  role  of  numerical  simulations  in  these  boundary  and  device  studies  has 
been  to  act  as  surrogates  for  measurements  that  are  not  feasible.  In  one  case, 
obtaining  cathode  boundary  fields  from  measurements  was  not  possible.  Thus 
for  long  devices  the  sensitivitiy  of  the  numerical  results  to  numerical  changes  in 
the  boundary  conditions,  when  coupled  to  experiments,  provided  the  key  to  the 
role  of  contacts  on  device  behavior  (1).  For  submicron  devices,  the  difficulties 
of  direct  correlation  of  experiment  with  specific  transport  phenomena  are 
apparent  and  simulation  through  parametric  studies  will  provide  a  key  to  the 
role  of  boundaries.  But  the  description  of  transport  on  a  submicron  scale  is  still 
inadequate  and  the  descriptive  role  of  boundaries  is  correspondingly  weak.  For 
example,  most  space  charge  dependent  problems  still  treat  the  background  as  a 
“jellium”  distribution.  The  discrete  nature  of  impurities  is  ignored,  as  are 
structural  variations  in  the  contacts.  The  extent  to  which  this  affects  such 
measurements  as  current-versus-voltage  is  yet  to  be  determined.  Notwithstand¬ 
ing  these  uncertainties,  a  considerable  amount  of  information  can  be  obtained 
by  extrapolation  from  the  ideal  cases  which  can  provide  bounds  on  the  limits 
of  transport  through  both  the  boundary  and  active  region  of  the  device. 
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ABSTRACT 

Recent  results  concerning  spatial  and  temporal  transport  in 
subraicron  devices  identify  significant  aspects  of  the  role  of 
boundary  conditions,  scaling  for  suggesting  new  materials,  and  struc¬ 
tural  device  changes.  These  results  are  discussed  as  a  means  of 
achieving  high  speed  and  high  frequency  devices. 


INTRODUCTION 

There  are  several  recent  results  concerning  spatial  and  temporal 
transport  in  submicron  devices  that  are  likely  to  have  an  impact 
on  the  design  of  future  high  frequency  sources.  These  results,  which 
emerge  from  Monte  Carlo,  and  momentum  moment  equation  solutions  to  the 
"Wigner-Boltzmann"  quantum  transport  equation^ 
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and  the  Boltzmann  transport  equation  (BTE) 
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identify  crucial  aspects  of  boundary  conditions,  the  role  of  scaling 

in  choosing  suitable  materials,  and  the  significance  of  alterations 

in  otherwise  simple  device  structures  for  achieving  high  speeds. 

These  tonics  are  reviewed  below.  (It  is  noted:  In  equation  (1)  the 

position  gradient  in  the  brackets  operates  only  on  the  potential 

energy,  f  is  a  single  coordinate  and  momentum  distribution  function 
w 

f  (x,p)  =  -  ^  f dyT*  (x  +  T  (x  -  (3) 

w  ZTin  f  l  z 

and  Y(x)  represents  the  state  of  the  system  in  the  coordinate  repre¬ 
sentation)  . 


SCATTERING  MODIFICATIONS  TO  BALLISTIC  TRANSPORT 

Ballistic  transport  implies  carrier  transport  unimpeded  by  inter 
actions  (electrostatic,  or  otherwise)  with  other  carriers,  or  with 
scattering  events.  The  extent  to  which  scattering  centers  are  sensed 
by  transiting  electrons  is  therefore  of  significance.  The  first  set 
of  calculation  shown.  Figure  1  (Ref.  2),  represents  scattering  events 
in  GaAs  for  a  collection  of  electrons  entering  a  uniform  field  region 
with  an  initial  energy  of  approximately  0.30ev.  Intervalley  P-L  ener 
separation  for  this  Monte  Carlo  calculation  is  0.33ev.  It  is  seen 
that,  with  the  exception  of  very  high  fields,  approximately  50%  of  th 
carriers  are  unscattered  over  the  first  500  A. 

The  velocity  versus  distance  curves  for  this  calculation  are  dis 
played  in  Figure  2  (Ref  2),  and  it  is  seen  that  high  speeds  over  use¬ 
ful  distances  can  be  achieved  even  in  regions  where  a  high  number  of 
scattering  events  has  occured.  The  optimum  conditions  for  this  ap¬ 
pear  to  require  moderate  fields,  generally  near  the  threshold  field 
for  electron  transfer,  and  moderate  injection  energies.  The  depend¬ 
ence  on  the  latter  is  displayed  in  Figure  3  (Ref. 2),  where  the  lowest 
achievable  velocities  over  a  distance  of  1500° A  are  for  entry  elec¬ 
trons  with  zero  initial  velocity.  It  is  also  seen  that  an  optimum 
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OlSTANCe  (A) 


Figure  1.  Percentage  of  unscattered  electrons  versus  distance. 

Entrance  electrons  have  a  finite  energy.  From  reference 
2,  with  permission. 

entry  energy  exists;  for  energies  near  the  F-L  separation,  scatter¬ 
ing  quickly  reduces  the  net  drift  velocity.  The  optimum  conditions 
identified  by  these  calculations  imply  that  voltage  control  near 
Che  entrance  boundary  must  be  near  30mv. 

The  examples  of  Figures  1  through  3  are  for  uniform  fields,  and 
carriers  subject  to  sudden  changes  in  field.  Studies  in  which  spatial 
gradients  accompany  nonuniform  fields  are  more  recent,  several  of  which 
are  considered  at  this  Workshop.  One  aspect  is  considered  below. 


SPATIAL  AND  TEMPORAL  TRANSIENTS 

A  relatively  direct  way  of  handling  spatial  and  temporal  trans¬ 
ients  is  through  moments  of  the  transport  equations  (1)  and  (2). 
Little  has  been  done  with  the  moments  to  the  Wigner  -  Boltzmann  equa¬ 
tion  and  the  following  deals  exclusively  with  BTE  moments.  The 
moments  of  the  BTE  form  an  infinite  hierarchy  of  moment  equations. 
Each  equation  introduces  a  higher  order  moment  not  defined  by  the 
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Figure  2.  Velocity  versus  distance  for  electrons  with  a  finite 

entrance  energy  subject  to  a  sudden  change  in  electric 
field.  From  reference  2,  with  permission. 


given  set  of  balance  equations.  The  most  common  form  of  the  balance 
equations  is  olit.iitu’d  by  assuming  a  distribution  function  of  the 
form 


f  -  f  .  _  i  exp  -  {  m(u-v)  )  ..  . 

,,,3/2  2  k  T.  ^ 

I  .  O  1 

I 

for  each  species.  More  generally  (Ref.  A) 

^2  3 

^  \rnir,  9u,‘^\lm  Du.Du  (5) 

k  k  1  Kim 

where  it  is  insisted  that  the  first  terra  in  the  expansion  gives  the 
correct  local  values  of  density,  velocity  and  energy.  The  coeffic¬ 
ients  in  the  above  expansion  are  inversely  proportional  to  density 
and  are  model  dependent.  The  first  correction  to  the  local  equilib¬ 
rium  balance  e<|uation.s,  the  so-called  hydrodynamic  approximation  yielc 
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OlSTANCe.  (A) 


Figure  3.  Velocity  versus  distance  for  electrons  with  varying  en¬ 
trance  energies  subject  to  a  sudden  change  in  electric 
field.  From  reference  2,  with  permission. 
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where  =  n^m^v.,  ,  o  is  a  stress  tensor 

arising^from^nonuniform  velocity  distributions,  and  k  is  the  thermal 
conductivity.  The  stress  term  is  dissipative  in  that  when  a  nonuni¬ 
form  velocity  distribution  is  impressed  on  an  electron  stream  there 
will  be  reactive  forces  tending  to  smooth  them  out.  Detailed  informa¬ 
tion  regarding  suitable  values  for  these  terms  are  not  available,  and 
they  are  regarded  as  phenomenological  entries;  they  and  their  modifi¬ 
cation  by  the  collision  terms  are  to  be  studied. 


div  (<  grad  T) 


PHYSICS  AND  MODELING  CONSIDERATIONS  FOR  VLSI  DEVICES 


69 


The  boxed  terms  above  represent  contributions  arising  from  the 
nonspherical  nature  of  the  distribution  function.  The  underlined 
terms  are  ignored  in  the  drift  and  diffusion  approximation.  Figure 
4  shows  the  result  of  a  solution  to  the  balance  equations  for  two 
level  transfer  for  a  two-terminal  device  whose  structure  is  receiving 
considerable  attention  as  an  illustration  of  velocity  overshoot.  The 
structure  considered  here  consists  of  a  low  doped  region  L'v^lym 
sandwiched  between  two  regions  of  hloher  donine  concentration.  The 
structure  is  subject  to  a  bias  of  1  volt.  Four  interfaces  are  involved 
but  most  attention  is  focused  on  the  two  interior  faces  where  signifi¬ 
cant  gradients  in  field  occur,  all  associated  with  differences  be¬ 
tween  the  background  and  mobile  carrier  density.  The  carrier  con¬ 
centration  is  displayed  in  figure  4c,  for  the  case  of  injecting  cathode 
contacts  (Ref.  4).  It  is  seen  that  (a)  significant  charge  injection 
is  present  in  the  low  doped  region,  and  (b)  electron  transfer  occurs 
near  the  downs tream^sect ion  of  the  low  doped  region  and  is  significant 
in  the  downstream  N  region.  It  is  noted  that  the  field  dependence  of 
carrier  density  in  the  T  valley,  Np(F)  is  not  the  same  as  that  asso¬ 
ciated  with  uniform  field  steady  state  values.  Spatial  gradients  re¬ 
sult  in  a  spatial  lag. 

Also  shown  in  figure  4c  is  the  drift  velocity  of  carriers  in  the 
r-valley,  and  it  is  apparent  that  the  carriers  acquire  speeds  consid¬ 
erably  in  excess  of  the  steady  state  peak  velocity  of  electrons  in 
GaAs.  These  velocities  are  not,  however,  much  different  in  value  than 
those  associated  with  F  valley  electrons  under  uniform  field  condi¬ 
tions.  The  velocity  shown  in  figure  4c  is  not,  however  the  same  as 
the  carrier  velocity  computed  under  uniform  field  conditions.  The 
velocities  computed  here  include  temperature  as  well  as  momentum 
gradients,  and  are  properly  defined  as  the  current  flux/carrier 
density.  Thus,  increased  velocity  here  is  often  associated  with  de¬ 
creased  carrier  density,  through  current  continuity.  The  key  to  all 
device  design  is  that  high  velocities  must  be  accompanied  by  only  mar¬ 
ginal  transfer  out  of  the  high  mobility  section  of  the  conduction 
band.  This  will  require  device  lengths  shorter  than  the  active  region 
of  the  device  shown  in  this  figure.  (Modifications  to  this  statement 
arising  from  variations  in  material  parameters  are  considered  in 
Ref.  4). 

The  calculations  of  figure  4  show  the  situation  where  a  high 
carrier  velocity  occurs  near  the  downstream  edge  of  the  low  doped 
region.  Often  high  velocities  are  required  near  the  upstream 
boundary  as  indicated  by  figures  1  through  3.  Figure  5  (Ref.  6) 
shows  a  spatially  dependent  result  for  carriers  subjected  to  a  pre¬ 
specified  value  of  electric  field.  The  calculation  is  for  silicon  and 
the  highest  carrier  velocities  occur  near  the  0.4vjm  boundary.  An 
accompanying  plot  of  carrier  density  (Ref.  6,  figure  9)  shows  the 
lowest  level  of  carrier  density  in  the  region  of  highest  velocity. 
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Figure  AC.  Fractional  population  and  central  value  velocity.  (Fn 
Ref.  4.,  with  permission.) 


Coupling  the  above  results  to  the  need  for  moderate  fields  for 
high  speed  operation  a  useful  multiterminal  device  very  likely  will 
require  the  presence  of  a  local  moderate  field  region  that  "kicks" 
electrons  into  high  speed  regions.  The  simplest  high  field  in¬ 
jecting  region  could  come  from  local  charge  depletion  or  "notches" 
(Fig.  6).  Note  that  in  Figure  6  the  gate_^is  treated  generically. 

It  could  be  a  metal  Schottky  contact,  a  p  AlGaAs  layer,  etc. 


SCALING 

In  the  absence  of,  or  in  concert  with  structural  variations  in 
devices  it  is  also  necessary  to  examine  changes  in  material  paramet* 
The  change  in  material  parameters  should  have  as  its  goal  high  mobij 
ity  and  high  characteristic  velocities.  Thornber  (Ref.  8)  provided 
general  set  of  guidelines  for  choosing  material  parameters  through 
an  alteration  of  the  scattering  rates,  one  of  which  is  duscussed  be! 


The  collision  term  in  the  Boltzmann  transport  equation  is 
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Figure  6.  Schematic  of  a  high  energy  injection  region  in  a  gated 
submicron  structure. 
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The  one  illustrated  below  involves  uniformly  altering  the  scattering 
rate  by  a  constant  F,  and  results  in  the  following  relation 


v^(x,t,F)  =  v(rx,rt,  F/r)  (ii) 

Where  the  right-hand  side  is  the  unsealed  velocity.  In  this  scaling 
mobility  is  altered  but  saturated  drift  velocity  is  not. 


A  dramatic  consequence  of  the  effects  of  scaling  is  illustrated 
below,  but  it  is  necessary  to  note  that  the  possibility  of  relevant 
material  scaling  over  a  broad  range  of  electric  field  values  is  re- 
mote.  Rather,  it  is  more  likely  to  be  appliable  over  a  restricted 
range  of  field  values.  This  is  illustrated  for  GaAs.  For  GaAs,  over 
a  field  range  of  approximately  3  ->■  15  kv/cm,  intervalley  F-L  coupling 
is  a  dominant  scattering  mechanism.  Figure  7  shows  the  velocity  field 
characteristic  for  F-L-X  ordering  with  three  different  values  of  the 
F-L  deformation  potential.  All  these  valleys  are  taken  as  parabolic 
with  while  all  other  material  parameters  are  those  of  Little¬ 

john,  et  al.  (Ref. 9).  (We  note  that  similar  values  for  material  con¬ 
stants  yield  differences  between  Monte  Carlo  and  balance  equation  cal¬ 
culations.)  In  Figure  7  the  fields  at  points  1,  2  and  3  occur  at  valu 
5.6  kv/cm,  8.2  kv/cm  and  12  kv/cm,  respectively.  The  dominant 
F  valley  toial  momentum  scattering  rates  are  4 . 78x10^ ^/sec . ,  7.32x10^^ 
sec.,  and  10.0x10^ ^/sec. ,  respectively.  The  ratios  F/F  are 
1.17x10  ^^kv-s/cm,  1.12x10  ^^kv-s/cm  and  1.20x10  ^^kv-s/cm,  respect¬ 
ively  in  general  agreement  with  the  rule  of  Equation  11. 


Figure  8  illustrates  the  consequences  of  uniform  scaling  on  veloc 
ity  overshoot.  The  solid  curve  displays  ovei.hoot  for  a  gallium 
arsenide  element  subject  to  a  field  of  27kv/cm  (Ref.  10).  The  dashed 
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Figure  7.  Steady  state  velocity  versus  field  for  T-L-X  ordering  as 

obtained  from  balance  equations.  Different  curves  are  for 
different  values  of  the  deformation  potential  coupling 
constant  for  intervalley  F-L  scattering. 


(dotted)  curve  is  a  sketch  of  a  scaled  curve  for  r-2(r=0.5)  and  a 
field  of  S'ikv/cm  (13.5kv/cm).  For  r=2,  high  overshoot  velocities 
occur  over  a  short  period  of  time  and  require  high  bias  levels.  By 
reducing  T  to  values  below  1.0  high  overshoot  velocities  are  retained 
for  a  longer  period  of  time  and  at  lower  values  of  field  than  for  GaAs. 
This  result  is  highly  significant  insofar  as  it  suggests  higher  transit 
velocities  for  extended  temporal  scales,  and  indicates  a  direction  of 
material  selection  for  high  frequency  sources. 


GENERAL  CONCLUSION 

It  has  been  known  since  1969  that  the  compound  semiconductors, 

GaAs  in  particular,  are  theoretically  capable  of  providing  high  fre¬ 
quency,  near  lOOGHZ  oscillations.  Overshoot  has  also  been  known  since 
that  time  (Ref.  11).  The  difficulty  in  attaining  high  frequency  three- 
terminal  and  sometimes  two-terminal  operation  has  over  the  past  decade 
been  attributed  to  inadequate  contact  regions  (Ref.  12),  material 
preparation,  etc.  Many  of  the  high  frequency  problems  were  thought 
to  be  reduced  by  going  to  small  dimensions,  where  in  addition  to 
shorter  transit  lengths  the  benefits  of  overshoot  would  emerge.  While 
it  is  still  too  early  to  state  how  these  benefits  can  be  implemented 
in  practice,  it  is  clear  that  special  contacts  or  injection  regions 


74 


H.  L.  GFHUBIN  ET  AL. 


Figure  8.  Transient  velocity  characteristics  for  scaled  and  unsealed 
scattering  rates.  Solid  curve  is  from  reference  10. 


are  required.  It  appears,  therefore,  that  the  1969  questions  as  to 
why  devices  are  not  operating  closer  to  theoretical  limits  is  still 
valid  today.  The  dual  solution  approaches  of  the  last  decade  also 
appear  equally  valid  today,  where  the  pursuit  of  novel  device  struc¬ 
tures  emphasizing  the  boundary  role  is  continuing,  while  simultaneous 
new  material  directions  are  sought. 
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In  this  letter,  the  relationship  between  the  characteristic  function  for  two  arbitrary  noncommuting  observables  and  a 
generalized  Wigner  distribution  function  is  established.  Ihis  distribution  function  is  shown  to  have  no  simple  interpretation 
in  the  sense  of  probability  theory  but,  in  lieu  of  its  special  properties,  can  be  used  directly  for  calculating  the  expectation 
values  of  observables. 


Whereas  classical  transport  physics  is  based  on  the 
concept  of  a  probability  distribution  function  which 
is  defined  over  the  phase  space  of  the  system,  in  the 
quantum  formulation  of  transport  physics,  the  con¬ 
cept  of  a  phase  space  distribution  function  is  not  pos¬ 
sible  inasmuch  as  the  noncommutation  of  the  posi¬ 
tion  and  momentum  operators  (the  Heisenberg  uncer¬ 
tainty  principle)  precludes  the  precise  specification  of 
a  point  in  phase  space.  However,  within  the  matrix 
formulation  of  quantum  mechanics,  it  is  possible  to 
construct  a  “probability”  density  matrix  which  is  of¬ 
ten  interpreted  as  the  analog  of  the  classical  distribu¬ 
tion  function. 

There  is  yet  another  approach  to  the  formulation 
of  quantum  transport,  based  on  the  construction  of 
the  Wigner  distribution  function  [1].  As  we  shall 
show,  this  distribution  function  has  no  simple  inier- 
pretation  in  the  sense  of  probability  theory  [2]  but, 
in  lieu  of  its  special  properties,  can  be  used  directly 
for  calculating  expectation  values  [3-5]  of  observ- 

'  Partially  supported  by  the  Army  Research  Office. 


ables  in  a  manner  quite  analogous  to  that  of  classical 
theory,  i.e.  by  integrating  the  product  of  the  observ¬ 
able  and  the  Wigner  distribution  function  over  all 
phase  space. 

This  letter,  in  part,  reviews  the  salient  features  of 
the  Wigner  distribution  function.  Although  the  Wigner 
function  is  generally  defined  in  terms  of  all  the  gener- 
aUzed  coordinates  and  momenta  of  the  system  in 
question  as 

P^^x^...x„,P^...p„)  =  ~  f  dy,...d>-„ 

X +5  3'l.-->r„  (>) 

X  'I'(x,  _  ...jr„  -?>•„) 

X  exp|i(p,y,  +  ...  * P„y„)/fl]  , 

we  will  discuss  the  properties  of  the  Wigner  function 
in  terms  of  a  single  coordinate  and  momentum.  In 
this  case,  we  let 


Volume  87A,  number  4 


PHYSICS  LtTTtRS 


4  January  1982 


where  'l'(jc)  refers  to  the  state  of  the  system  in  the  co¬ 
ordinate  representation. 

The  distribution  function  of  eq.  (2)  has  interesting 
properties  in  that  the  integration  of  this  function  over 
all  momenta  leads  to  the  probability  density  in  real 
space;  conversely,  the  integration  of  this  function  over 
all  coordinates  leads  to  the  probability  density  in  mo¬ 
mentum  space.  In  mathematical  terms, 


function  to  the  classical  distribution  function  is  appar¬ 
ent  by  examining  the  equation  of  time  evolution  for 
Pv/^x,p).  Upon  assuming  that  'l'(x)  in  eq.  (2)  satisfies 
the  Schiodinger  equation  for  a  system  with  hamilto- 
nianW  =  p2/2m  +  K(x),  it  can  be  readily  shown  that 
P}f/(x,p)  satisfies  the  equation 

a/’w/a/  +  (p/m)a/>^/9x-te-/’^=o,  (S) 

where 


S  (-iy> 

n  =  0 


(2/1  1)1 


/  =  {3a) 

and 

f  P\\/ix,p)<ix -<l)'{p)<li{p) ,  (3b) 

where 

=  J  e-‘P^/«T'(A:)dx  . 

It  follows  immediately  from  eq.  (3)  that,  for  an 
observable  W[x,p)  which  is  either  a  function  of  mo¬ 
mentum  operator  alone  or  of  position  operator  alone, 
or  any  additive  combination  therein,  the  expectation 
value  of  the  observable  is  given  by 


It  is  evident  that  in  the  limit  Pi-*  0,9  •  in  eq.  (6) 
becomes 

e-p^  =  ~idyidx)(dp^ldp)  (7) 

so  that  eq.  (5)  reduces  to  the  classical  collisionless 
Boltzmann  equation. 

The  Wigner  distribution  function  defined  in  eq.  (2) 
is  derivable  [6)  from  the  Fourier  inversion  of  the  ex¬ 
pectation  value  (with  respect  to  state  'P(x)  of  the  op¬ 
erator  (here,.*  and  p  satisfy  the  commuta¬ 

tion  relation  [x,p]  =  ih).  As  such, 

P)  =  -~  If  0)e-HTP-9*dT  dd  ,  (8a) 

4»i'‘ 

where 


<lv^>  =  JJ WPy)^(x,p)  d*  dp  , 


(4)  Cw(t,0)= /'l'*(*)ci(^P-"»^>'l'(*)d*  ,  (8b) 


which  is  analogous  to  the  classical  expression  for  the 
average  value.  Herein  lies  the  interesting  aspect  of  the 
Wigner  distribution  function;  the  result  of  eq.  (4) 
suggests  that  it  is  possible  to  transfer  many  of  the  re¬ 
sults  of  classical  transport  theory  into  quantum  trans¬ 
port  theory  by  simply  replacing  the  classical  distribu¬ 
tion  function  by  the  Wigner  distribution  function. 
However,  unlike  the  density  matrix,  the  Wigner  dis¬ 
tribution  function  itself  cannot  be  viewed  as  the  quan 
turn  analog  of  the  classical  distribution  function  since 
it  is  generally  not  positive  definite  and  noiiunique 
[/’^y(*,p)  of  eq.  (2)  is  not  the  only  bilinear  expres¬ 
sion  (1 ,3-5]  in  '1'  that  satisfies  eq.  (3)]. 

Further  resemblance  of  the  Wigner  distribution 


and  the  interval  of  integration  is  (-'*>,  “>)  unless  other¬ 
wise  specified.  In  order  to  sliow  that  the  right-hand 
side  of  eq.  (8a)  is  indeed  the  Wigner  distribution  func¬ 
tion  as  defined  in  eq.  (2),  note,  from  the  Baker- 
Hausdorff  theorem  |7j,  that  ed/P'*®*)  can  be  written  a 

pi(T/7+9i)  =  giT;7/2giOigiTp/2 

in  wluch  case  Cyv(t,  0)  of  eq.  (8b)  becomes 


Cw(t.0)=  r  le-'^P/2q/(x)j*e'»^[e-‘^P/2vK*)]  d*  , 

(10) 


which  further  reduces  to 
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Cw(t,0)=  f  -\Th)e'^^<l'(x+{Th)dx  .  (11) 

Then,  by  inserting  C\v(t,  8)  of  eq.  (1 1)  into  the  right- 
hand  side  of  eq.(8a),  integrating  over  the  variable  6 
by  using  the  relation 

J  eie(x'-x")  ^2it6(x'  -  x”) , 

and  letting  t  =  the  desired  result  is  obtained. 

The  method  outlined  above  to  arrive  at  the  Wigner 
distribution  function  is  based  on  the  notion  of  a  char¬ 
acteristic  function.  The  characteristic  function  of  an 
observable,  A ,  with  respect  to  state  !'}')  (here,  the 
Dirac  notation  is  utilized  for  purposes  of  generality) 
is  defined  as 

C4(?)=  <'I'|e‘Wlvl/>  ,  (12) 

where  ^  is  a  real  parameter.  Assuming  >1  to  possess  an 
eigenvalue  spectrum  given  by  A  \A  ')  =A'\A'},  C^(U 
can  be  evaluated  in  the  /T-representation  as 

C^(?)=  fdA'f  dA"  ('\'\A')<A'  \A")<A"  m  . 

Since  (A'\c'^'^  \A")  =  e'^-^  5(A'  -  A")  in  the/l'-rep- 
rcsentation,  (|)  in  eq.  (13)  reduces  to 

=  ,  (14) 

where  I'l',^'!'  =  IC-l'l'P)!'  ^  1'{A‘),  the  probability 
distribution  function  for  measuring  .-1 '  while  in  slate 
I'D.  Hence,  the  characteristic  function  lot  A  is  the 
Fourier  transform  of  the  probability  disinbiiiion 
function  P(A').  Subsequent  inversion  of  eq.  ( 14) 
above  leads  to 

P(A‘)  =  ~-  d%.  (1.S) 

The  Wigner  distribution  function  was  derived  by 
taking  the  Fourier  transform  of  the  characteristic 
function  for  e'^rp*ex)  yjg^  yg  (l,g  connection  be¬ 
tween  the  probability  distribution  function  and  the 
characteristic  function  for  a  given  observable,  this  ap¬ 
proach  seems  to  be  a  natural  way  of  obtaining  a  dis¬ 
tribution  function  for  momentum  and  position.  Un¬ 
fortunately,  the  noncommutative  nature  of  the  two 


observables  destroys  the  convenient  probability  inter¬ 
pretation  of  the  characteristic  function  implicit  in 
eq.(15). 

In  order  to  demonstrate  this  point,  assume  the  char¬ 
acteristic  function  for  two  noncommuting  observables, 
A  and  B,  to  be 

C/1B«1«2)=  -  (16) 

Observables  y4  and  B  are  assumed  to  have  eigenvalue 
spectra 

A[A')=A'IA')  ,  B\B')  =  B'\B’)  ,  (17) 

and  are  chosen  so  that  [A,  [A,  B]  ]  =  [B,\A,  fi]  ]  =0. 
This  assumption  is  imposed  so  that  the  identity 

ei(tizl+{2B)  e'^2®e-fiOl-4.fll/2  (jg) 

may  be  used. 

Inserting  eq.  (18)  into  eq.  (16)  while  obtaining  the 
matrix  elements  of  e‘UA  jn  the  /I '-representation  and 
e*^J®  in  the  ^'-representation  results  in 


X  e‘Ui''l'-"Ej^')<'F|/l'>U'|5'>(5'|'I'>  . 

In  eq.  (19),  it  is  assumed  that  \A ,  B]  is  a  c-number  in¬ 
dependent  of  the  eigenvalues /I  ’  and  B'.  We  define 
F(A\B'),  the  generalized  Wigner  distribution  func¬ 
tion  * ' ,  to  be 

FiA‘,B')={'\>\A'){A'\B')iB'\'^)  ,  (20) 

so  that 

F(/1',B')=J -  fd^t  Tdb  cUUM‘.fil/2 

Ctr)^  ■ 

(21) 

It  is  evident  from  c<)s.  (20.21 )  that 


The  form  of  the  generalized  distribution  funclion  derived 
IS  sensitive  to  the  manner  in  which  cxp|i(t)/4  *  tjB))  is 
expanded.  I  or  example,  if  instead  of  eq.  (18),  we  used  the 
foim5exp(it2fi)cxp(iti/4)exp(t,tj|/l,  ej/2),  exp(i{,>t/2) 
X  expdtjfllexpdt  ,/i/2).  or  expdtjB/Z)  expdt)/4)  exptitj 
X  BID,  all  of  which  are  equivalent,  we  would  indeed  ob¬ 
tain  a  different  form  of  the  generalized  distribution  func¬ 
tion.  yet  one  which  obeys  the  sum  rules  of  eq.  (22). 
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and 

/ F(A’.  B')  dB' 


(22a) 


(22b) 


Thus,  eq,  (21 )  establishes  the  relationship  between  the 
characteristic  function  for  two  arbitrary  noncommut- 
ing  observables  and  the  generalized  Wigner  distribution 
function.  The  generalized  distribution  function  has 
the  essential  properties  of  the  conventional  Wigner 
function  in  that  an  integration  of  the  generalized 
function  over  the  eigenvalue  spectrum  of  one  observ¬ 
able  leads  to  the  probability  density  in  the  canonically 
conjugate  observable  [eq.  (22)]. 

There  is  no  simple  probability  interpretation  of 
F{A',  B')  in  eqs.  (20,  2 1 )  because  of  the  necessary 
overlap  between  the  states  of  the  noncommuting  ob¬ 
servables.  However,  if  A  and  B  are  made  to  commute 
so  that  M')  and  16')  are  a  common  set  of  eigenvec¬ 
tors,  then  F{A' ,  B')  reduces  to  the  probability  distri¬ 
bution  function  for  A  and  B. 

Finally,  it  is  noted  that  the  conventional  Wigner 
distribution  function  for  observables  A  and  B  is  ob¬ 
tained  from 

(2^)“ 

X  (23) 

with  $2)  defined  in  eq.  (16),  whereas  the  alter¬ 

native  distribution  function,  F(A',  B').  introduced  in 


eqs.  (20,  21)  differs  from  the  Wigner  fuiiction  due  to 
the  presence  of  the  phase  factor  e^'  ^2 1/2  in  the 
integrand  of  eq.  (21).  For  >4  =x  and  6  =  p,P\i/{x,p) 
in  eq.  (23)  reduces  to  the  Wigner  function  of  eq.  (2), 
whereas  F(x,p)  defined  from  eq.  (20)  becomes 

F(x,p)  =  ^/  dy>i'-(x)^ix-y)e'Pylf> 

where  0(p)  is  defined  in  eq.  (3b).  It  is  evident  that 
there  is  a  family  of  functions  which  are  bilinear  in  iF 
yet  satisfy  the  sum  rules  of  eqs.  (3a,  b). 

There  are  some  interesting  questions  to  be  resolved 
concerning  the  uniqueness  and  positive  definiteness 
of  Wigner-type  quantum  distribution  functions.  Never¬ 
theless,  these  distribution  functions  serve  a  useful 
purpose  for  calculating  quantum  mechanical  observ¬ 
ables  in  transport  [5]  studies  and  numerous  solid- 
state  [8,9]  problems. 
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ABSTRACT 

The  size  of  exi^ected  future  very-small  devices  will  result  in 
their  being  strongly  coupled  to  the  environment  in  which  they  are 
located.  In  this  pajjer,  we  examine  several  of  the  limitations  on  the 
device  physics  that  can  be  expected  to  arise  in  these  structures, 
including  the  role  of  the  environment  and  its  effect  on  ballistic 
transport.  In  addition,  we  look  at  submicron  arrays  of  devices  in 
a  lateral  surface  superlat t ice  (LSSL)  and  examine  transport  in  such 
an  LSSL  tlirough  a  Wignei'  function  and  Monte  Carlo  approach. 

INTRODUCTION 

Many  people  have  examined  the  limit  to  which  semiconductor  de¬ 
vices  can  be  scaled  downward*”^.  While  small  devices  in  the  range 
0.1  -  0.3  pm  gate  length  have  been  made**,  problems  such  as  intercon¬ 
nections,  electromigration  and  thermomigration  of  metallization,  and 
power  density  within  the  device  strongly  affect  the  packing  density 
and  device  size  that  can  be  achieved^.  Yet,  no  one  has  evaluated  the 
operation  and  performance  of  very-small  semiconductor  devices,  i.e. 
those  that  can  be  conceived  in  the  sub-0.1  pm  size  range,  and  the 
interactions  within  arrays  of  such  devices.  One  reason  for  this  is 
that  the  transport  within  such  a  device  cannot  be  treated  in  isolation. 
Because  of  the  size  of  such  a  very-small  device,  it  is  coupled  strongly 
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to  the  environment  in  which  it  is  located.  The  basic  transport  equa¬ 
tions  cannot  be  separated  from  their  casual  boundary  conditions  and 
both  of  these  factors  must  be  modified  to  account  for  the  influence 
of  the  device  environment^.  In  previous  work,  some  of  us  have  pre¬ 
sented  a  formalism  to  address  this"**^,  but  this  formalism  remains 
untried  due  to  its  inherent  complexity.  In  the  present  work,  we  wish 
to  examine  qualitatively  several  of  the  limitations  on  the  device 
physics  that  can  be  expected  to  arise  in  the  very-small  semiconductor 
device,  particularly  with  respect  to  the  finiteness  of  source  and 
sink  regions,  contact  regions  and  barriers,  and  the  influence  of  these 
latter  quantities  on  the  ballistic  transport  that  may  exist  in  the 
device®’^.  In  particular,  we  examine  these  effects  through  the  use 
of  discrete  area-preserving  phase-space  maps  which  display  the  char¬ 
acteristics  of  transport  in  a  generic  device.  In  addition,  we  treat 
the  cooperative  interactions  that  can  arise  in  densely  packed  arrays 
of  devices. 


BALLISTIC  TRANSPORT 

Potential  barriers  within  the  device  can  play  a  significant  role 
in  the  quantum  ballistic  transport  of  carriers  through  it.  Such 
barriers  are  found  in  very-small  devices,  for  example,  to  confine 
carriers  to  the  active  region®,  and  are  an  intimate  part  of  devices 
such  as  the  planar-doped  barrier  transistor^®,  tunneling  barrier  de- 
vices^^,  real-space  transfer  devices^®,  or  superlattice  avalanche 
photodiodes^®.  When  a  barrier  is  present  at  the  contact  region,  care 
must  be  taken  to  adequately  handle  turning-point  reflection  of  the 
electrons  from  this  barrier.  Even  when  the  electron  has  sufficient 
energy  to  pass  over  the  barrier,  there  is  a  well-known  quantum  mechan¬ 
ical  reflection  at  the  barrier  interface.  If  the  potential  barrier  is 
smooth,  i.e.  introduces  a  transition  over  many  wavelengths,  the  re¬ 
flection  and  wave  function  matching  can  be  handled  by  well-known 
approximation  techniques  such  as  the  WKB  approximation''*,  in  which 
the  potential  barriers  represent  turning  points  for  a  near-classical 
path.  In  the  very-small  device,  however,  the  barriers  are  expected 
to  be  sharp  on  the  scale  of  the  electron  wavelength  and  care  must  be 
exercised  in  matching  wavef unct ions  and  determining  reflections.  The 
reflection  problem  is  further  complicated  in  real-space  transfer  de¬ 
vices'®  due  to  the  different  band  structure  on  either  side  of  the 
barrier.  Here,  additional  terms  arise  due  to  the  spatial  variation 
of  the  effective  mass. 

If  the  potential  barriers  are  slowly  varying  on  the  scale  of  the 
wave  packet,  the  trajectories  are  largely  those  of  the  classical 
motion.  Even  if  this  is  not  the  case,  as  we  expect  for  the  very-small 
device,  nearly  semiclassical  trajectories  can  be  expected  if  the  var¬ 
iation  of  the  action  is  limited  to  a  few  low-order  derivatives'®. 

In  semi-classical  systems,  the  phase  space  of  the  classical  motion 
forms  a  natural  framework  in  which  to  examine  problems  such  as  these. 
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While  classical  one-dimensional  transport  appears  to  be  basically 
simple,  there  exists  recent  work  that  suggests  that  even  this  simple 
problem  contains  a  number  of  unexpected  subtleties In  this  paper, 
we  begin  to  Investigate  the  role  of  barriers  in  submicron  devices 
through  the  use  of  such  finite,  area-preserving  phase-space  mappings 
of  classical  dynamics.  The  use  of  such  mappings  reveals  a  variety 
of  complicated  structure^^* and  has  had  recent  success  in  explain¬ 
ing  the  cause  of  excess  noise  in  Josephson  junction  parametric  amp¬ 
lifiers^®*  The  results  that  we  obtain  indicate  that  if  these 

mappings  are  applicable  to  the  VLSI  scale,  then  present  concepts  of 
submicron  transport  may  require  substantial  generalization. 

The  types  of  sublet'es  to  which  we  are  referring  are  best  ill¬ 
ustrated  by  Fig.  1.  There,  we  are  using  an  area-preserving  mapping 
of  particle  position  (horizontal)  and  momentum  (vertical)  within  a 
device  active  region  bounded  by  two  Gaussiam  potentials  V(x)  (over¬ 
laid  in  the  figure).  In  addition,  an  electric  field  has  been  applied 


Figure  1.  The  phase-spacing  mapping  of  an  initial  uniform  distri¬ 
bution  of  electrons  in  a  generic  very-small  semiconductor 
device.  The  contact  regions  are  represented  by  a  twin- 
gaussian  potential  V(x).  The  horizontal  direction  is 
position  while  the  vertical  direction  is  momentum. 
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The  canonical  mappings  are  area-preserving  since  we  are  considering 
a  collisionless  (conservative)  system  and  looking  at  the  ballistic 
transport.  The  classical  differential  equations  of  motion  are 


3V(x) 

P  =  -  ‘ 

X  =  +  p/m  , 

and  the  discreet  classical  area-preserving  maps  are 


(1) 


3V 

Pn+1 

=  Pn 

- 

X  .  1 

n+l  , 

X  .  , 

=  X 

+  T  0  /m  - 

n+l 

n 

n 

In  Fig.  1, 

the 

potential  has  been  scaled  so  that  the  two 

(2) 


have  a  weak  overlap.  For  small  values  of  the  total  energy,  resonant 
orbits  occupy  the  central  part  of  the  figure.  These  orbits  are  in  the 
region  of  the  classically  integrable  motion.  For  larger  values  of 
the  energy,  however,  the  orbits  are  such  that  the  particles  are  swept 
out  of  the  well  by  the  field.  An  energy  dissipating  collision  can 
drop  a  high  energy  particle  into  the  resonant  region,  thus  trapping 
it  within  the  structure.  We  expect  these  particles  to  contribute  a 
diffusive  component  of  current.  Large  angle  scattering,  however,  can 
move  a  particle  to  a  back-flowing  orbit,  which  effectively  causes  a 
reflection  of  particles  from  the  device  input. 


A  number  of  interesting  factors  arise.  The  potential  is  generic, 
in  the  sense  that  it  is  similar  in  form  not  only  to  the  devices  men¬ 
tioned  above,  but  also  to  the  sinusoidal  force  term  in  the  Josephson 
junction  devices^®’ In  the  latter  case,  interaction  between  devices 
can  lead  to  a  parametric  pumping  of  the  potential  which  yields  period¬ 
doubling  bifurcations  and  chaotic  behavior.  The  system  (2)  is  not 
truly  generic  though,  as  the  set  (x,p)  are  not  proper  action-angle 
variables.  We  have  examined  the  behavior  through  potential  pumping, 
with  dissipation,  through  replacing  the  first  equation  of  (2)  with 


dv 


n+1. 


-  T 


9x 


Typ  +  TF  sin(Qt  ) 
no  n 


(3) 


where  Vq  is  the  set  of  Gaussian  potentials.  The  factor  Y  is  an  ef¬ 
fective  damping  factor.  In  Fig.  2,  we  plot  the  (Fo,fi)  plane  results. 
The  curve  is  a  separatrix  below  which  a  stable  device  results.  Above 
the  curve,  the  device  is  unstable.  No  period-doubling  bifurcations 
within  a  single  device  are  found,  contrary  to^®. 
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Figure  2.  The  separatrix  in  (Fq,Q)  space  for  parametric  pumping  of 
the  device  of  Fig.  1.  For  values  of  (F^.f^)  above  the 
curve,  electrons  are  ejected  from  the  device,  while  for 
values  below  the  curve  the  device  is  stable. 


DEVICE  ARRAYS 

When  device  sizes  begin  to  shrink  toward  the  0.1pm  or  less  region, 
the  line-to-line  capacitance  in  dense  device  arrays  begins  to  dominate 
the  total  node  capacitance^.  This  parasitic  capacitance  leads  to  a 
direct  device-device  interaction  outside  of  the  normal  circuit  or 
architectural  design.  In  conventional  descriptions  of  LSI  circuits, 
each  device  is  assumed  to  behave  in  the  same  manner  within  the  total 
system  as  it  does  when  it  is  isolated.  In  the  dense  arrays  discussed 
here,  this  will  no  longer  be  the  case. 

The  possible  device-device  coupling  mechanisms  are  numerous  and 
include  such  effects  as  the  capacitive  coupling  mentioned  above,  but 
also  include  such  effects  as  wave  function  penetration  or  tunneling 
and  charge  spill-over.  Formally,  however,  one  may  describe  these 
effects  on  system  and  device  behavior  by  assuming  the  simplest  form 
of  coupling.  Arrays  of  devices,  interacting  in  this  manner,  form  a 
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lateral  surface  superlattice^^ * ^ ^ .  Lateral  superlattices,  in  wh 
the  superstructure  lies  in  a  surface  or  heterostructure  layer,  o 
considerable  advantages  for  obtaining  superlattice  effects  in  pi 
technology.  While  a  surface  MOS  structure  is  formally  similar  t- 
array  of  CCD  devices^® » ,  superlattices  can  also  be  fabricated 
the  use  of  electron  and  ion  beam  lithography  and  selective  area  < 
taxial  growth^ If  the  coupling  is  capacitive,  then  the  limita 
to  a  spacing  less  than  the  de  Broglie  wavelength  is  removed**.  Wi 
have  examined  transport  in  such  lateral  surface  superlattices  (L! 
through  a  Wigner  function  approach.  Before  proceeding,  however, 
is  worth  noting  that  the  circuit  theory  view  of  LSSLs**  is  generii 
that  of  cellular  automata^ Many  of  the  image  processing  applic 
tions  proposed  for  LSSLs^**  arise  from  the  "games"  aspect  of  celli 
automata^^. 

Superlattice  structures  give  rise  to  sinusoidal  energy  mini¬ 
bands  with  relatively  narrow  widths.  The  shape  of  such  bands  re* 
in  interesting  electrical  transport  properties.  Lateral  surface' 
superlattices  have  cosinusoidal  bands  in  two  dimensions,  with  the 
third  dimension  quantized  with  discrete  energy  levels. 

We  have  calculated  the  transport  properties  for  such  LSSLs  i 
a  Wigner  function  formulation.  A  complete  integral  equation  is  c 
tained  for  the  Wigner  function  and  must  be  solved  to  obtain  the  t 
port  coefficients^^ However,  the  form  of  the  solutions  can  b 
found  by  taking  the  average  velocity  and  energy  from  the  first  an 
second  moments  of  an  equivalent  Wigner  representation  of  the  trar 
port  equation^®.  A  constant  electric  field  is  applied  in  the  pla 
of  the  sinusoidal  bands.  For  simplicity,  a  constant  relaxation  t 
T  can  be  used.  Solving  these  two  equations  simultaneously  leads 
the  velocity  and  energy  as  functions  of  the  field^® 

(27t/fiL)^eF  T 

<v^>  -  (e+<Ej_>  -  <E^>)  1  +  (2iTeF  t/fiL)^  ’ 

z 

<E  >  -  e  -  <E  > 

<E>  =  e  +  <Ej^>  + 

z 

where  L  =  2it/D,  D  is  the  LSSL  spacing,  e  is  one-half  the  band-wid 
and  <Eq>  and  <E]^>  are  the  equilibrium  and  transverse  energies,  re 
spectively . 

The  velocity  as  a  function  of  field  has  the  same  basic  analy 
expression  as  that  obtained  by  Lebwohl  and  Tsu®®  except  for  the  e: 
factor  In  front  of  the  expression  for  the  field.  The  difference 
in  effect  caused  by  the  different  equilibrium  distribution  chosen 
In  the  latter  paper,  the  authors  assumed  the  initial  equilibrium 
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tribution  is  very  close  to  a  zero  temperature  Fermi-Dirac  function, 
while  here  the  distribution  is  a  finite  temperature  Wigner  distribution 
and  Includes  the  physics  of  the  energy  band  shape.  Note  that  the 
velocity  as  a  function  of  the  electric  field  shows  the  negative 
differential  mobility  predicted  earlier  and  also  exhibits  the  general 
shape  expected  of  the  velocity-field  curve. 

To  more  exactly  illustrate  the  mobility  and  negative  differential 
conductivity,  we  have  carried  out  a  Monte  Carlo  calculation.  The 

O  1 

scattering  processes  are  calculated  for  the  model  of  lafrate  et  al.  , 
and  thus  are  for  the  system  of  GaAs/GaAlAs.  The  scattering  rates  for 
acoustic  and  polar  optical  phonons  have  been  obtained  using  a  two 
dimensional  donsity-of-states  for  cosinusoidal  energy  bands.  A  Van 
Hove  sing,ulnrity  occurring  in  the  density-of-states  produces  a  sing¬ 
ularity  in  the  scattering  rate,  which  was  removed  by  including  the 
sel  f-ener;'y  corrections  due  to  the  phonons  in  the  vicinity  of  the 
singularity.  The  widely-spaced  discrete  energy  levels  in  the  third 
dimension  allows  scattering  and  transport  in  that  direction  to  be 
nc)’,  lec  toil . 

In  this  surface  superlattice,  as  in  others,  the  conduction  band 
splits  into  subbands.  Here  the  lowest  energy  subband  was  nearly  flat. 
Therefori',  transport  dominantly  occurs  in  the  next  higher  minibands. 

The  satellite  valleys  and  next  subband  are  at  energies  of  0.2eV  and 
0.3eV,  respectively,  above  the  subband  considered.  Their  contribution 
to  the  transport  of  electrons  is  insignificant  since  there  are  no 
intermeil late  energies  through  which  the  electrons  can  scatter  to  aid 
pc^pulation  of  upi)er  bands. 

Tlie  overall  transport  properties  of  this  system  are  calculated 
by  .an  ensemble  Monte  Carlo  technique.  The  results  of  the  simulation 
.are  the  ve  1  oe  i t y-f  ield  curves  and  are  shown  in  Fig.  3.  The  lower 
curve  lesults  for  a  field  applied  along  one  of  the  (10)  basis  vectors 
of  the  stpiare  lattice  array  of  cylinders  while  in  the  top  curve  the 
fielti  i  .apjilied  along  a  (11)  direction.  At  low  fields,  both  curves 
show  .1  liiu’.ar  region  as  expected  for  most  structures.  At  approxi- 
m.ately  .S-10  kV/cm  the  curves  begin  to  bend  over  to  the  peak  near 
13  kV/cm.  As  the  field  is  further  increased,  the  velocity  begins  to 
decrease  .and  for  this  model  continues  to  decrease  to  zero  as  the 
fielii  tends  to  infinity. 


In  summary,  transport  and  scattering  in  a  generic  surface  super¬ 
lattice  structure  exhibits  a  negative  differential  mobility  arising 
from  Bloch  oscillations^ ^ .  This  surface  superlattice  negative  dif¬ 
ferential  mobility  is  expected  to  be  useful  at  much  higher  frequencies 
than  that  due  to  the  conventional  Gunn  effect.  Alternatively,  related 
instabilities  may  be  an  ultimate  limit  on  very  large  scale  integra¬ 
tion  ^  . 
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FIELD  (kV/cm) 


Figure  3.  The  drift  velocity  (as  a  function  of  electric  field)  for 
a  LSSL  at  300  K.  The  GaAs/GaAlAs  model  of  ref.  21  has 
been  used. 
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